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Abstract 

The  control  of  the  transition  from  laminar  to  turbulent  flow  in  a  boimdary  layer 
of  a  flat  plate  is  investigated  using  numerical  simulations.  The  numerical  model  is 
based  on  the  incompressible  Navier-Stokes  equations,  which  are  coupled  with  the 
energy  equation  through  the  temperature  dependent  viscosity.  A  fully  implicit  finite- 
difference  spectral  method  was  used  to  solve  the  governing  equations.  The  numerical 
model  allows  for  the  spatial  evolution  of  the  disturbances  in  a  non-paraUel  boundary 
layer.  Active  control  of  wave  packet  disturbances  in  the  non-isothermal  boundary 
layer  is  studied  in  detail  Wave  packet  disturbances  are  created  in  the  flow  field  by 
simulating  the  effect  of  thermally  activated  heater  elements  on  the  plate  surface. 
Through  a  controlled  spanwise  variation  of  the  temperature  of  the  heater  elements, 
two-  and  three-dimensional  wave  packet  disturbances  can  be  studied.  The  propaga¬ 
tion  and  amplification  of  the  wave  packet  disturbances  in  the  boundary  layer  is 
examined.  The  heater  elements  on  the  plate  surface  act  as  locally  strong  heat  sources 
causing  thermal  wakes  within  the  boundary  layer  that  spread  in  the  downstream 
direction.  A  transfer  function  technique  is  used  for  the  control  strategy.  The  transfer 
function  is  based  on  the  vortidty  response  to  a  finite  temperature  fluctuation  at  the 
heater  strip  and  is  obtained  from  the  numerical  simulations.  With  additional  heater 
segments  (controller)  located  downstream  of  an  excitation  source,  the  possibility  of 
attenuating  wave  packet  disturbances  is  investigated.  With  the  numerical  transfer 
function,  a  successful  control  strategy  for  the  wave  packet  cancellation  could  be 
devebped.  Initially,  for  the  low  amplitude,  two-dimensional  disturbances  in  the  tran¬ 
sition  process  the  wave  packet  disturbances  could  be  almost  completely  cancelled.  For 
the  attenuation  of  three-dimensional  wave  packet  disturbances,  the  transfer  function 
technique  was  extended  to  allow  for  spanwise  variations.  The  attenuation  of  three- 
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dimensional  wave  packets  with  the  modified  transfer  function  technique  was  almost 
equally  as  successful  as  for  the  purely  two-dimensional  flow  disturbances.  For  the 
simulation  of  the  three-dimensional  flow  development  with  no  control  applied,  non¬ 
linear  interaction  of  wave  components  of  the  wave  packet  first  appeared  for  the 
oblique  modes  in  the  low  frequency  range,  which  was  also  observed  in  experimental 
investigations.  The  attenuation  of  only  the  two-dimensional  components  of  a  three- 
dimensional  wave  packet  disturbance  delaj^  the  onset  of  the  nonlinear  interaction  of 
the  oblique  spanwise  modes  in  the  lower  frequency  range. 
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1.  Introduction 


The  study  of  laminar  to  turbulent  transition  for  wall  bounded  shear  flows  is  of 
great  practical  interest  for  the  application  of  flow  control.  The  delay  of  transition 
from  laminar  to  turbulent  flow  at  high  Reynolds  numbers  reduces  the  friction  drag 
on  streamlined  bodies.  The  effective  control  of  laminar  flow,  for  example  to  extend 
the  regime  of  laminar  flow  on  the  surface  of  airplane  wings,  may  lead  to  an  increased 
range  and  speed,  or  reduced  fuel  costs  of  an  airplane.  On  the  other  hand,  in 
processes  where  the  mixing  of  the  fluid  projjerties  is  desirable  in  order  to  improve 
momentum,  energy  or  mass  transport,  an  acceleration  of  transition  to  the  turbulent 
state  is  advantageous.  Applications  of  the  control  of  laminar  flow  were  reviewed  in 
great  detail  by  Bushnell  and  Malik  (1985),  and  more  recently,  by  Wagner  et  al.  (1989) 
The  present  numerical  investigation  employs  the  flat  plate  model  of  a  boundary  layer 
flow  in  water.  The  instability  mechanisms  leading  to  turbulence  are  reviewed  first  to 
gain  insight  into  the  relevant  issues  prior  to  the  discussion  of  laminar  flow  control. 


One  route  of  transition  to  turbulence  begins  with  the  occurrence  of  small  ampli¬ 
tude,  two-dimensional  disturbance  waves  (Tollmien  Schhchting  waves).  These  distur¬ 
bance  waves  propagate  in  the  direction  of  the  flow.  The  development  of  these  distur¬ 
bances,  that  is  amplification  or  decay,  can  be  closely  predicted  by  the  linear  stability 
theory.  At  later  stages  the  Tollmien  Schlichting  (TS)  waves  develop  a  spanwise 
periodic,  three-dimensional  deformation,  and  associated  with  this  is  the  formation  of 


longitudinal  vortices.  In  the  experiments  of  Klebanoff  et  al.  (1962),  this  behavior  was 


enhanced  by  placing  physical  geometric  disturbances  on  the  plate  surface  at  selected 


spanwise  intervals.  A  very  regular  peak-valley  structxire  appeared  (regions  of 


I 


17 


enhanced  and  reduced  amplitudes)  with  the  peaks  corresponding  to  regions  of  max¬ 
imum  wave  amplification  and  the  valleys  to  regions  of  minimum  amplification. 

The  growth  rates  at  the  peak  positions  in  the  Klebanoff  pattern  are  much  larger 
than  the  original  TS  growth  rates  and  subsequently  a  localized  high  shear  layer  is 
formed.  Then,  shortly  after,  highly  nonlinear  processes  lead  to  breakdown  and  to 
turbulent  spots,  which  spread  to  form  a  fully  turbulent  boundary  layer. 

Unlike  the  above  mentioned  route  to  transition,  at  somewhat  lower  amplitude 
levels  (between  0.2%  and  1%  of  the  free  stream  velocity),  Saric  and  Thomas  (1983), 
and  Kachanov  and  Levchenko  (1984)  observed  a  different  pattern.  Instead  of  an 
aligned  pattern  seen  in  the  experiments  of  Klebanoff  et  al.  where  peaks  followed 
peaks,  and  valleys  followed  valleys,  this  alternative  pattern  consisted  of  a  staggered 
structure,  i.e.  peaks  followed  valleys  and  vice  versa. 

The  staggered  structure  at  the  lower  initial  forcing  amplitudes  leads  to  a  strong 
amplification  of  the  subharmonic  mode,  and  also  to  the  growth  of  a  broad  band  of 
lower  frequency  disturbances.  Kachanov  and  Levchenko  (1984)  did  not  observe  the 
characteristic  high  frequency  stages  of  the  Klebanoff  type  transition.  In  addition  to 
the  above  mentioned  routes  to  tiubulence,  certain  bypass  mechanisms  may  occur, 
that  is,  if  the  external  forcing  results  in  large  enough  amplitudes  in  the  initial  stages, 
nonlinear  effects  may  trigger  breakdown  to  turbulent  flow  in  a  more  direct  manner. 

In  an  early  theoretical  study  on  the  breakdown  process,  Craik  (1971)  proposed  a 
weakly  nonlinear  resonance  model  based  on  the  interaction  of  a  two-dimensional  and 
two  three-dimensional  oblique  modes  with  opposite  signs.  The  resonance  interaction 
would  involve  modes  whose  phase  velocities  match  in  streamwise  direction.  In  this 
theory,  which  is  based  on  the  temporal  amplification,  this  would  typically  involve 
three-dimensional  wavenumbers  =  1/2  Oio-  However,  the  relatively  large  distur- 
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bances  observed  in  Klebanoffs  experiment  could  not  be  explained  by  Grade's  theory. 

Much  later  the  experiments  by  Saric  and  Thomas  (1983),  and  Kachanov  and 
Levchenko  (1984)  provided  clear  evidence  of  the  existence  of  the  subharmonic  mode. 
In  the  smoke  wire  flow  vistializations  by  Saric  and  Thomas  two  different  peak-valley 
structiu'es  were  identified.  For  a  very  low  initial  disturbance  amplitude  (  .3%  based 
on  the  free  stream  velocity)  the  spanwise  wavenumbers  matched  Craik's  theory  well 
A  spanwise  to  streamwise  wavenumber  ratio  of  y  /  a  =  .67  was  reported.  However,  at 
slightly  higher  initial  forcing  amplitudes  (.4%),  a  spanwise  to  streamwise 
wavenumber  ratio  of  y/ a  =  1.46  was  identified,  which  cannot  be  explained  with 
Craik's  theory. 

An  explanation  of  this  difference  was  later  given  by  Herbert  (1984).  He  studied 
the  stability  of  the  Blasius  profile  and  a  superimposed  plane,  time-periodic  T5  wave 
with  respect  to  three-dimensional  Squire  modes.  According  to  Herbert's  calculation 
the  growth  of  spanwise  wave  components  can  be  explained  by  the  fact  that  the  phase 
velocities  of  the  Squire  modes  (disturbance  modes  with  zero  normal  disturbance 
amplitude)  are  independent  of  the  spanwise  wave  angle.  This  growth  of  spanwise 
modes  is  an  important  feature  of  the  boundary  layer.  The  effective  control  of  any 
three-dimensional  waves  in  the  boundary  layer  will  have  to  address  these  important 
aspects. 

However  for  important  practical  applications,  such  as  flow  over  airplane  wings, 
the  instability  mechanism  that  ultimately  leads  to  transition  to  turbulent  flow  is 
essentially  linear  at  first,  before  nonlinear  interactions  take  over.  In  addition,  for  a 
disturbaxKe  background  of  small  amplitudes  this  linear  regime  covers  by  far  the  larg¬ 
est  distance  in  the  entire  transition  process.  Therefore  the  first  control  attempt  of 
transition  will  have  to  concentrate  on  the  early  transition  stages. 
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1.2  Linear  Theory 

In  natural  transition  ToUmien  Schlichting  waves  can  be  observed  when  the 
amplitudes  of  the  background  disturbances  are  very  small.  As  long  as  the  amplitudes 
of  the  disturbances  in  the  boimdary  layer  remain  small  one  can  study  the  stability  of 
steady  laminar  boimdary  layers  by  decomposing  the  flow  quantities  as 
c{=  qo(x)  +  where  qo  is  the  base  flow,  and  q'ix4)  is  a  small  disturbance  quan¬ 

tity,  the  momentum  equations  can  be  linearized.  The  resulting  equations  can  be 
further  simplified  by  assiuning  a  localty  parallel  flow,  so  that  terms  involving  the 
mean  normal  velocity  and  streamwise  derivatives  of  mean  quantities  can  be  omitted. 
Then,  by  using  local  velocity  profiles,  the  stability  of  a  boimdary  layer  can  be  investi¬ 
gated  (cf.  Schlichting,  1982,  p454  ff). 

A  solution  ansatz  in  the  form  (J>  (y)  •  e  •<“*  -P*)  leads,  after  some  algebra,  to 
the  derivation  of  the  Orr-Sommerfeld  equation 

(U-c  ,  (1 .0) 

outeji 

where  a,  P  are  the  complex  wavenumber  and  the  frequency,  respectively,  and  (l)(y)  is 
the  complex  streamfunction,  which  is  a  function  only  of  the  distance  normal  to  the 
wall 

The  nondimensional  wavenumber  a  and  the  Reynolds  number  Rej^  are  based  on 
the  displacement  thickness,  and  the  prime  denotes  differentiation  with  respect  to  y. 
This  equation  together  with  the  boundary  conditions 

=  =  0  ;  y=0,oo  (1.1) 

forms  an  eigenvalue  problem  with  the  parameters  a,  c,  and  with  the  eigenfunction 
(y).  For  example,  a  real  valued  P  and  a  complex  a  =  a+io,  represents  the  spatial 
amplification  case,  where  -<x,>0  denotes  spectrally  amplified  and  -<x,<0  spatially 
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damped  waves.  For  each  a  fixed,  and  Uiy)  ,U  "(y)  given  from  the  Blasius  solu¬ 
tion,  the  eigenvalue  c  can  be  computed.  A  finite  number  of  distinct  eigenvalues  exists, 
together  with  a  continuous  spectrum  for  which  the  eigenfunctions  vary  sinusoidally 
as  (Jordinson,  1971,  and  Grosch  and  Salwen,  1978).  However  the  solution  of 

the  stability  problem  for  the  boundary  layer  yields  only  one  discrete  amplified  mode 
for  an  unstable  case.  All  the  other  discrete  modes  are  damped. 

The  neutral  stability  curve,  where  ot,-  =  P,-  =  0,  for  a  zero  pressure  gradient 
boimdary  layer  is  shown  in  Figure  3.  One  can  discern  an  upper  and  a  lower  branch 
of  the  neutral  stability  curve,  which  is  divided  by  the  point  where  the  tangent  to  the 
curve  is  normal  to  the  Re^-axis.  This  point  defines  a  critical  Reynolds  number  Recrit  > 
because  for  lower  Res^  no  amplified  oscillations  can  exist.  This  minimum  Re^rit  for 

the  parallel  flow  assumption  is  approximately  520.  Between  the  lower  and  upper 
branches  of  the  neutral  curve,  the  disturbances  amplify  exponentially,  and  outside  of 
these  branches  the  disturbances  decay. 

A  great  deal  of  theoretical  and  experimental  effort  has  been  expended  on  the 
study  of  single  frequency  disturbances.  However,  natural  excitation  by  free  stream 
disturbances  produces  modulated  wave-trains  or  wave  packets,  which  are  observed 
to  break  down  to  turbulence  in  a  much  more  dramatic  manner  than  two-dimensional 
single  fi'equency  waves. 
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1.3  Wave  Packets 

Often  natural  transition  from  laminar  to  turbulent  flows  in  a  boundary  layer  is 
initiated  by  broad  band,  pulse-Kke  disturbances  or  wave  packets.  Such  wave  packets 
can  already  be  seen  as  modulated  waves  in  the  time  traces  at  certain  streamwise  sta¬ 
tions  in  the  early  experiments  of  Schubauer  and  Skramstad  (1948).  For  such  distur¬ 
bances  in  the  flat  plate  boundary  layer,  the  amplified  modes  vary  continuously  with 
respect  to  one  or  more  frequencies,  and  the  most  amplified  mode  is  merely  one 
among  equals  because  the  neighboring  modes  have  relative  growth  rates  and  phase 
velocities  close  to  the  most  amplified  one.  The  amplified  modes  represent  a  group  of 
waves,  with  their  development  being  in  time  and  space. 

Gaster  and  Grant  (1975),  in  an  experiment  on  the  development  of  isolated  wave 
packets,  kept  the  initial  amplitudes  of  a  localized  wave  packet  small  to  study  the 
linear  evolution  of  the  packet,  and  to  avoid  nonlinear  distortion  of  the  boundary 
layer  in  the  early  transition  stages.  Thus,  the  evolution  of  the  packet  is  controlled  by 
the  amplification  rate  and  the  dispersion  relation  of  constituent  waves.  The  wave 
packet  first  followed  the  linear  prediction  manifested  in  the  Gaussian  envelope  of  the 
amplitudes.  At  locations  further  downstream,  these  envelopes  became  distorted, 
which  may  indicate  a  non-uniformity  of  the  mean  flow,  or,  the  appearance  of  non¬ 
linear  interactions  of  the  wave  packet  modes. 

Gaster  (1975)  then  applied  a  theoretical  model  based  on  the  solution  of  the  Orr- 
Sommerfeld  equation.  By  a  direct  summation  of  the  eigenfunctions,  he  tried  to 
assemble  the  wave  packet  disturbances  of  the  experiment.  The  initial  conditions  for 
the  summation  of  the  eigenfunctions  were  chosen  to  be  close  to  that  of  the  experi¬ 
ment.  The  agreement  between  experiment  and  theory  was  remarkably  good,  except 
for  the  development  of  a  two  peak  phenomena  (a  distortion  of  the  envelope)  that 
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was  observed  in  the  experiment  at  locations  further  downstream  from  the  excitation 
source. 

In  a  direct  numerical  simulation  with  the  complete  Navier-Stokes  equations, 
Fasel  (1983)  has  calculated  two-dimensional  wave  packets  with  an  initially  large 
amplitude  input,  to  demonstrate  that  the  distortions  of  the  envelopes  in  the  experi¬ 
ments  of  Caster  and  Grant  (1975)  maybe  attributed  to  nonlinear  effects.  These  simu¬ 
lations  showed  a  broadening  of  the  amplitude  spectrum  due  to  wave  interactions  of 
the  amplified  modes  in  the  wave  packet. 

Recently,  Konzelmann  (1991)  simulated  numerically  the  experiment  of  Caster 
and  Grant  (1975)  using  the  complete  three-dimensional  Navier-Stokes  equations.  The 
parameter  of  this  computation  closely  matched  those  of  the  experiment  and  a 
remarkably  good  agreement  between  the  numerical  simulation  and  the  experiment 
was  fotmd.  These  simulations  provided  strong  evidence  that  the  breakdown  process 
was  subharmonic.  This,  for  example,  would  explain  the  appearance  of  the  low  fre¬ 
quency  components  for  the  oblique  modes. 

1.4  Control 

Hrst  attempts  to  control  the  laminar  boimdary  layer  concentrated  on  passive 
methods.  There,  the  modification  of  the  mean  velocity  profile,  for  example,  using 
pressure  gradients,  wall  suction,  constant  heating  or  cooling  of  the  wall,  alters  the 
stability  characteristics  of  the  base  flow.  For  the  flat  plate  boundary  layer  in  water, 
the  eff^  of  constant  cooling  is  destabilizing,  whereas  the  effect  of  a  constant  heating 
is  stabilizing.  Since  the  viscosity  of  water  at  ambient  temperature  decreases  sharply 
with  increase  in  temperature,  heating  yields  a  fuller  velocity  profile  which,  according 
to  the  linear  stability  theory,  results  in  a  more  stable  flow.  In  the  water  experiments 
of  Strazisar  et  al.  (1977),  uniformly  heating  the  boimdary  layer  increased  the 
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mininrnm  critical  Reynolds  number  and  shifted  the  maximum  amplification  rates  to 
lower  frequencies.  Strazisar  et.  aL  tested  several  wall  temperatures  in  their  experiment 
and  showed  that  the  transition  to  turbulence  could  be  delayed  considerably. 

Lowell  (1974) ,  and  Wazzan,  Okamura  and  Smith  (1968)  extended  the  linear  sta¬ 
bility  theory  (see  eqn.  1.0)  to  the  non-isothermal  case.  They  allowed  for  a  variation  of 
fluid  viscosity  with  temperature,  and  Lowell  also  introduced  an  additional  stability 
equation  for  the  temperature.  Both  linear  models  predicted  the  disturbance  behavior 
observed  in  the  experiments  only  qualitatively. 

In  a  numerical  simulation,  Bestek  et  al.  (1987)  used  the  full  Navier-Stokes  equa¬ 
tions  coupled  with  an  energy  equation  and  showed  qualitatively  good  agreement  to 
the  experimental  results  of  Strazisar  (1968).  Their  numerical  model  allowed  the  study 
of  a  non-uniform  as  well  as  a  vmiform  temperature  distribution  along  the  plate  sur¬ 
face  and  showed  the  stabilizing  as  well  as  the  destabilizing  effect  of  heating  and  cool¬ 
ing,  respectively,  of  the  boimdary  layer  in  water. 

Economic  considerations  limit  the  practical  value  of  passive  heating  control.  The 
large  amount  of  heating  necessary  to  reduce  the  fluid  viscosity  in  water  affects  a 
technical  operation  economically  because  of  increased  operation  and  maintenance 
cost  for  the  heating  devices.  Thus,  more  recently,  the  application  of  control  of  lam¬ 
inar  flow  has  focussed  on  active  methods.  For  this  control  approach  the  disturbance 
waves  are  directly  attacked,  their  amplitude  attenuated  or  enhanced  (depending  on 
the  application),  while  the  base  flow  stability  characteristics  are  approximately  unal¬ 
tered.  Thus,  the  amplification  rates  remain  imchanged,  but  the  disturbance  waves 
begin  to  grow  again  from  a  much  lower  amplitude  level.  The  basic  principle  of 
active  control  is  the  superposition  of  two  disturbance  waves.  With  proper  phase  and 
amplitude  control  between  disturbance  waves  of  the  same  frequency,  reinforcement 
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or  cancellation  of  the  disturbance  wave  can  be  observed. 

In  experimental  investigations  several  techniques  were  applied  for  an  active  con¬ 
trol  such  as:  heated  wires  (Corke,  1989),  heater  strips  (Liepmann  et  al.,  1982a,b), 
vibrating  ribbons  (Thonuis  1983  ,  Pupator  and  Saric  1989,  and  Milling  1981)  and  flexi¬ 
ble  walls  (Wehrmann,1965).  For  example,  with  the  heater  strips  used  in  the  experi¬ 
ments  by  Liepmann  et.  al.,  controlled  TS  waves  could  be  produced.  In  these  experi¬ 
ment  it  was  shown  that  active  control  of  ToUmien  Schlichting  waves  in  a  water  boun¬ 
dary  layer  with  one  activating  and  one  controlling  heater  strip,  led  to  a  considerable 
increase  of  the  transition  Reynolds  number  in  comparison  with  the  uncontrolled  case. 
Also  Liepmann  et  al.  were  able  to  excite  wave  packets  using  the  heater  strip  tech¬ 
nique.  However,  their  experiment  exclusively  addressed  the  control  aspect  of  two- 
dimensional  single  frequency  disturbances. 

Active  control  of  two-dimensional  wave-trains  was  studied  experimentally  by 
Milling  (1981),  Thomas  (1983),  and  Pupator  and  Saric  (1989).  In  all  experiments  the 
disturbances  were  introduced  with  a  vibrating  ribbon  and  subsequently  controlled  by 
additional  ribbons  placed  inside  the  bormdary  layer  and  downstream  of  the  activator 
ribbon.  Milling  (1981)  and  Thomas  (1983)  concentrated  on  the  control  of  single  two- 
dimensional  disturbances,  and  Pupator  and  Saric  (1989)  focussed  on  the  control  of 
random  two-dimensional  disturbances.  In  all  experiments  the  amplitudes  of  the  two- 
dimensional  disturbances  were  considerably  reduced,  although  a  completely  distur¬ 
bance  free  environment  could  not  be  achieved.  Eventually  the  three-dimensional  dis¬ 
turbance  background  in  the  flow  caused  the  breakdown  to  turbulent  flow.  With  the 
low  free  stream  turbulence  wind  tunnel  of  Pupator  and  Saric,  the  natural  transition 
number  of  Re=3.5'10*.  With  the  control  of  the  two-dimensional  random  disturbances 
the  transition  Re-  number  increased  to  5  -10^.  Apparently,  the  control  of  two- 
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dimensional  disturbances  appreciably  increased  the  regime  of  laminar  flow,  yet  for 
further  improvement  a  two-dimensional  control  is  not  sufficient,  because  of  the 
three-dimensional  disturbance  amplification. 

Numerical  studies  on  active  transition  control  have  also  been  performed  by  Lau- 
rien  and  Kleiser  (1989)  and  McMurray  et.  al.  (1983).  They  considered  temporally 
growing  TS  waves  only,  and  were  limited  in  their  investigation  to  a  parallel  base 
flow.  In  both  investigations  a  wall  normal  disturbance  component  was  applied  on 
the  entire  plate  surface.  Using  the  wave  superposition  principle  the  disturbances  were 
attenuated  by  adjusting  the  amplitudes  and  the  phases.  However,  only  partial 
attenuation  of  the  two-dimensional  wave  was  achieved  and  it  was  concluded  that 
control  would  only  be  effective  if  applied  at  an  early  two-dimensional  stage. 

Recently  three-dimensional  spatially  amplified  disturbances  in  the  flat  plate 
boimdary  model  were  studied  numerically  by  Fasel,  Rist  and  Konzelmaim  (1987). 
Both  two-dimensional  and  three-dimensional  disturbances  were  introduced  into  an 
isothermal  boundary  layer  by  means  of  suction  and  blowing  through  a  slot  on  the 
plate.  The  motivation  for  these  simulations  was  to  investigate  the  mechanisms  of 
fundamental  and  subharmonic  disturbance  growth.  Fasel  et  al.  (1987)  achieved 
remarkably  good  agreement  with  the  experimental  investigations  of  Klebanoff  et  al. 
(1962)  for  the  fundamental  case,  and  of  Kachanov  and  Levchenko  (1984)  for  the 
subharmonic  case. 

Krai  (1988)  extended  the  numerical  approach  of  Fasel  et.  al  (1987)  to  model  dis¬ 
turbance  excitations  with  the  heater  strip  technique.  Krai  (1988)  studied  the  cancella¬ 
tion  of  two-  and  three-dimensional  spatially  amplified  disturbances  on  the  flat  plate 
boundary  layer.  However,  only  single  frequency  disturbances  were  considered  and 
emphasis  was  placed  on  the  control  of  the  subharmonic  and  fundamental  breakdown 
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processes.  Also  for  passive  control,  the  effect  of  uniform  and  non-uniform  plate  tem¬ 
perature  on  the  growth  rates  of  the  disturbances  was  investigated.  The  passive  con¬ 
trol  method  effectively  altered  the  growth  rate  of  the  two-dimensional  mode  which 
then  delays  the  growth  of  the  three-dimensional  components. 

1.5  Problem  Statement 

Active  control  of  the  laminar  boundary  layer  with  heater  strips  is  investigated 
using  numerical  simulations  based  on  the  complete  Navier-Stokes  equations.  From 
experimental  observations  it  is  evident,  that  disturbance  waves  in  the  boundary  layer 
evolve  spatialty.  Spatially  propagating  disturbance  waves  can  be  created  with  a  con¬ 
trolled  temperature  variation  over  a  heater  strip  as  was  demonstrated  experimentally 
by  Liepmann  et  aL  (1982a).  Such  temperature  variations  over  a  heater  strip  result  in 
a  local  change  of  the  fluid  viscosity,  and,  when  the  temperature  varies  with  time,  pro¬ 
duces  a  vortidty  fluctuation  on  the  wall  surface  and  its  immediate  vicinity.  In  a 
botmdary  layer  flow  of  water  along  a  flat  plate,  a  temperature  variation  of  20°C  at 
ambient  temperature  creates  a  vortidty  response,  which  corresponds  to  a  typical  dis¬ 
turbance  present  in  a  natural  transition  process. 

The  heater  strip  technique  also  has  the  advantage  of  not  placing  a  physical 
obstruction  in  the  flow  field,  imlike  vibrating  ribbons  or  heated  wires  placed  away 
from  the  wall,  inside  the  boundary  layer.  The  heater  strips  are  easily  activated  by  an 
alternating  current  which  produces  temperature  fluctuations  of  twice  the  current  fre¬ 
quency. 

For  the  present  investigation  a  numerical  method  is  used  that  is  based  on  the 
vortidty-  velodty  formulation  of  tiie  complete  Navier-Stokes  equations,  coupled  with 
the  energy  equation  through  a  temperature  dep)endent  viscosity.  The  simulation  are 
based  on  the  so-called  sp>atial  model,  thus  allows  for  the  spatial  amplification  of 
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disturbance  waves  in  a  laminar  boimdary  layer.  Finite  difference  approximations  are 
used  in  streamwise  and  wall  normal  directions/  and  spectral  approximations  are  used 
in  spanwise  direction  for  the  three-dimensional  simulations. 

Brst,  the  propagation  and  amplification  of  two-dimensional  wave  packet  distur¬ 
bances  are  investigated,  where  the  wave  packets  are  created  with  the  heater  elements. 
From  experiments  we  know  that  in  natural  transition  a  broad  band  of  frequency 
components  are  amplified.  Thus  numerical  simulations  of  wave  packet  disturbances, 
which  also  contain  a  broad  band  of  frequency  components  are  more  relevant  for 
investigating  the  control  of  natural  transition. 

A  control  strategy  based  on  a  numerically  determined  transfer  function  is  intro¬ 
duced  to  achieve  the  attenuation  of  wave  packet  disturbances.  The  transfer  function 
is  based  on  the  wall  vortidty  response  of  the  boundary  layer  from  a  pulse  tempera¬ 
ture  fluctuation  over  the  heater  strip.  The  applicability  of  the  control  strategy  is  first 
demonstrated  for  purely  two-dimensiortal  wave  packet  disturbances.  Then  the  control 
strategy  is  extented  for  attenuating  three-dimensional  wave  packet  disturbances.  The 
three-dimeitsional  wave  packet  distiubance  is  assumed  to  be  periodic  in  spanwise 
direction.  Then  the  wave  packet  can  be  spectrally  decomposed  with  respect  to  the 
spanwise  wavenumbers,  where  each  spanwise  mode  can  be  controlled  separately. 

For  a  practical  application  of  the  control  of  three-dimensional  disturbance 
waves,  the  heater  elements  have  to  be  divided  into  several  spanwise  segments,  which 
can  be  controlled  independently.  Such  arrays  of  spanwise  heater  elements  have  been 
successfully  used  in  experiments  by  Robey  (1986),  to  create  single,  oblique  distur¬ 
bance  waves  in  the  flat  plate  boundary  layer,  and  by  Nygaard  (1991),  to  force  three- 
dimensional  disturbances  in  the  botmdary  layer  of  a  splitter  plate,  where  the  distur¬ 
bances  then  propagated  into  a  free  shear  layer.  However,  both  experimental  invest!- 
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gations  did  not  attempt  the  cancellation  of  disturbance  waves.  ^ 

In  the  present  investigation  three-dimensional  control  is  also  attempted,  where 
spanwise  wavenumber  components  are  treated  individually.  Then  for  the  control  of 
the  individual  spanwise  wavenumber  components  different  control  procedures  can  be 
investigated.  For  example,  by  attenuating  only  the  dominant  two-dimensional  wave 
packet  and  leaving  the  spanwise  modes  xmcontroUed,  the  effect  on  the  three- 
dimensional  components  can  be  studied.  Also  the  control  of  wave  packet  distur¬ 
bances  for  the  two-  and  three-dimensional  wave  components  is  attempted.  The 
effectiveness  of  the  control  strategy,  applied  to  each  mode  individually,  is  examined. 
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2.  Governing  Equations 


The  flow  of  a  viscous  fluid  along  a  flat  plate  with  nonuniform  temperature  is 
governed  by  the  Navier-Stokes  equations  and  the  energy  equation.  In  order  to  study 
the  effect  of  a  temperature  field  on  the  flow,  the  Navier-Stokes  and  energy  equations 
must  be  coupled  through  temperature  dependent  fluid  properties.  The  present  inves¬ 
tigations  require  the  numerical  solution  of  the  governing  equations  for  both  two- 
dimensional  and  three-dimensional  disturbance  waves. 


Z1  Equations  for  Two-dimensional  Flow 

It  win  be  assumed  that  the  flow  is  incompressible,  the  fluid  is  Newtonian,  and 
that  the  spedfic  heat  capacity,  the  thermal  conductivity  and  the  density  are  indepen¬ 
dent  of  temperature.  These  assumptions  are  valid  for  water  and  many  other  liquids, 
as  long  as  in  the  investigated  control  volume,  large  pressure  and  temperature 
differences  can  be  avoided.  However  the  change  of  viscosity  with  respect  to  tempera¬ 
ture  cannot  be  assumed  constant  for  most  fluids.  For  water  at  20°C,  a  temperature 
increase  of  25®C  already  gives  rise  to  a  40%  change  in  the  kinematic  viscosity.  For 
the  derivation  of  the  Navier-Stokes  equations  it  will  be  assumed,  that  the  kinematic 
viscosity  is  a  function  of  temperature.  Since  the  temperature  is  expressed  as  a  scalar 
field  with  spatial  dimensions,  the  spatial  viscosity  gradients  in  the  Navier-Stokes 
equations  do  not  vanish. 

The  Navier-Stokes  equations  were  derived  with  the  nondimensional  variables 
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0  = 


T-T,  . 
T*-T«  ' 


(2.1) 


where  the  x,  y  coordinates  denote  the  streamwise  and  wall  normal  flow  directions 
(Figure  2),  respectively,  u,  v  are  the  streamwise  and  wall  normal  velocity  com¬ 
ponents,  (0  the  vortidty,  6  is  the  dimensionless  temperature,  u  the  kinematic  viscos¬ 
ity,  is  the  free  stream  velocity,  is  the  free  stream  temperature,  is  the  wall 
temperature,  L  is  a  characteristic  length,  and  u,,  is  the  free  stream  viscosity.  All  quan¬ 
tities  with  a  bar  are  dimensional. 


Then  the  Navier-Stokes  equations  in  dimensionless  form  can  be  derived  with  the 
above  nondimensional  variables  as; 
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(2.2b) 


The  viscosity  in  the  above  equations  (Z2)  is  assumed  to  be  a  function  of  the 
temperature  alone.  Consequently,  all  the  spatial  derivatives  which  involve  the  viscos¬ 
ity  could  be  replaced  with  the  chain  rule,  i.e.  3u/3x  =  (3u/30)  (30/3x ).  However  for 
clarity  in  the  derivation  of  the  vorticity  transport  equation,  the  spatied  viscosity 
derivatives  are  left  in  place. 


By  taking  the  curl  of  equation  Z2,  the  vorticity  transport  equation 
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can  be  derived,  where  the  vortidty  (o  is  defined  as 

By  Re  Bx 

To  solve  for  the  velocity  components,  two  Poisson  equations 

1  B^u  ^  3^  _  ^ 

Re  Bx^  By^  dy 

1  3^v  ^  3^  _  3cd 

Re  ax^  3y2  ~  3x 

are  used. 


(2.5a) 


(2.5b) 


The  energy  equation,  in  nondimensional  form  with  the  temperature  as  the 


dependent  variable,  is 


38  ...  38  ..  30  1  ,  1  3^0  .  320,  „ 

—  +  u  —  +  V  -r - ;r-(-;; - t  +  — 5")  =  0 

Bt  Bx  3y  Pr  Re  dx^  By^ 


The  viscous  dissipation  term  in  the  energy  equation  has  been  dropped,  since  the  Eck¬ 
ert  number  ,  Ec  =  /  (c  T  J,  for  incompressible  flow  at  ambient  temperature  is 

very  small. 

In  the  above  equations.  Re  =  (J^L/u,,  is  the  Reynolds  number,  Pr=x)jK^  is  the 
Prandtl  number,  where  u«  is  the  kinematic  viscosity  and  k„  the  thermal  diffusivity  at 
free  stream  temperature. 

For  the  kinematic  viscosity  an  empirical  relation  of  the  form 
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is  used,  where  Ta^=506  cuid  Tj,=  -140  are  empirical  constants  estimated  for  water 
(Truckenbrodt,  1981),  and  AT=T„-T^  is  the  difference  between  the  free  stream  tem¬ 
perature  and  a  wall  temperature. 

2.1.1  Boundary  Conditions 

In  this  work,  bovmdary  conditions  are  employed  that  Fasel  (1976)  had  success¬ 
fully  used  for  transition  simulations  for  the  flat  plate  boundary  layer.  Analogous 
conditions  are  applied  for  the  energy  equation.  In  the  following  discussion  of  the 
boimdary  conditions,  the  letters  A  to  D  refer  to  the  four  comers  of  the  integration 
domain,  which  is  shown  schematically  in  Figure  2. 

At  the  plate  surface,  A-B  in  Figure  2,  of  the  integration  domain,  the  conditions 
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are  iised,  where  the  interval  [X| ,  X2]  is  the  streamwise  extent  of  a  heater  strip.  Thus, 
no  slip  and  an  impermeable  wall  is  assumed.  The  relationship  (2.8d)  for  calculating 
the  wall  vorticity  was  derived  from  2.5b  using  the  continuity  equation  together  with 
the  no  slip  condition. 

Several  choices  of  boxmdary  conditions  to  model  the  heater  strip  are  available. 
Following  a  theoretical  analysis  for  thin-film  heated  elements,  which  can  be  used  to 
measure  the  local  wall  shear  stress,  an  approximation  for  the  spatial  temperature  dis¬ 
tribution  of  the  heater  strip  can  be  found.  Thin-film  heaters  are  based  on  the  Rey¬ 
nolds  analogy  of  heat-transfer  and  skin  friction.  The  heat  supplied  to  the  fluid  is 
related  to  the  local  value  of  the  skin  friction,  and  can  be  foimd  from  the  electrical 
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power  supplied  to  the  heated  film.  However,  one  part  of  the  electrical  power  sup- 
plied  to  the  heater  is  transferred  to  the  substrate  through  heat  losses.  In  a  theoretical 
analysis  of  thin-film  heating  elements.  Brown  (1967)  derived  an  expression  relating 
the  electrical  power  P  input  to  a  heated  film  and  the  wall  shear  stress  by 

=  ,  (Ma) 

with 


where  k  is  the  thermal  conductivity,  p  is  the  fluid  density,  p  is  the  dynamic  viscosity, 
W  and  L  are  the  spanwise  width,  and  the  streamwise  length  of  thin-film,  respectively, 
and  C  denotes  conductive  losses  to  the  substrata 


In  an  experiment  in  a  laminar  boimdary  layer  flow  along  a  flat  plate.  Brown 
(1967)  verified  the  relation  (2.9)  by  calibrating  a  heated  film  on  the  plate  surfoce.  The 
film  was  controlled  by  a  constant-temperature  feedback  bridge  and  values  of  P  were 
recorded  for  several  free  stream  velocities  and  wall  shear  stresses.  The  corresponding 
value  of  the  skin  friction  is  found  from  the  theoretical  solution  of  the  zero-pressure 
gradient  case.  The  experiments  supported  the  validity  of  equation  (2.9)  by  demon¬ 
strating  a  strong  linearity  between  the  electrical  power  input  and  the  cubic  root  of  the 
wall  shear  stress. 


The  derivation  of  equation  (2.9)  assumed  a  top-hat  temperature  distribution  for 
the  "physical"  streamwise  length  of  the  heated  film.  In  the  presence  of  heat  losses  to 
the  substrate  it  is  vinrealistic  to  assume  a  perfectly  insulated  film.  However,  Brown 
calculated  the  "effective"  streamwise  length  of  the  heated  film  from  the  calibration 
curve  by  measuring  the  slope  M  and  the  spanwise  width  W  and  using  equation 
(2.9b).  For  this  procedure  Brown  determined  the  effective  streamwise  length  of  the 
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heated  element  to  be  twice  the  actual  film  length. 

From  the  above  analysis  for  thin-film  heaters,  the  heater  strip  temperature  at  the 
plate  sur^e  in  the  interval  [xip:2]  is  approximated  to 

iclx—x  ) 

e^,(xj)  =  (1+T,(f ))  sin^(  ^  )  ,  (2.10) 

X2-X1 

where  (f )  is  a  heater  temperature  disturbance  signal. 

For  the  computation  of  the  steady  state  solution,  in  (ZIO)  is  set  to  zero.  The 
special  form  of  ,  which  introduces  the  wave  packet  disturbances  in  the  flow  field  is 
discussed  in  the  chapter  for  the  numerical  results. 

At  the  inflow  boundary,  A-D  in  Figure  2,  the  similarity  solution  (Blasius  solu¬ 
tion)  for  the  isothermal  boimdary  layer  is  imposed. 

U  —  ^Biatius  >  V  “  ^Blasius  >  (D  ®  ^Blasius  >  0  —  0  .  (2.11) 

The  flow  at  the  inflow  boimdary  is  assumed  steady  and  isothermal.  The  wave  packet 
disturbances  are  introduced  so  for  downstream  of  the  inflow  boimdary  that  no 
noticeable  upstream  effects  of  the  disturbances  are  present  near  the  inflow  boundary. 

At  the  top  of  the  integration  domain,  C-D  in  Figure  2,  the  flow  variables  are 
separated  to  a  steady  and  a  transient  part  according  to 

u=U  +  u'  ,  v  =  V  +  v',  (2.12) 


(0  =  £1  +  oa'  ,  0  =  0  +  0', 

where  U,  V,  Q,  0  are  the  steady  components,  and  v',  oa'  and  0'  denote  the  distur¬ 
bance  components.  The  following  conditions  are  used: 
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(0  =0  ,  0  =s0  . 

A  zero  pressure  gradient  in  streamwise  direction  is  assumed,  so  that  the  streamwise 
base  velocity  rennains  constant  at  the  upper  boundary  of  the  integration  domain.  The 
disturbance  components  are  assumed  to  decay  exponentially,  where  a  is  a  nondimen- 
sional  wavenumber  a  =  Q  L  of  a  most  amplified  frequency  component  of  the  wave 
packet. 

At  the  outflow  boimdary,  B-C  in  Figure  2,  the  conditions 
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are  used.  The  conditions  for  the  dishirbance  quantities  allow  a  wave  motion  in  the 
downstream  direction,  through  the  outflow  boundary  out  of  the  integration  domain. 
The  wavenumber  a  is  selected  from  the  maximum  amplified  disturbance  component 
of  the  wave  packet  at  the  outflow  boimdary. 


2.2  Equations  for  Three-Dimensional  Flow 

For  the  derivation  of  the  vortidty  transport  equation,  we  start  with  the  dimen¬ 
sional  Navier-Stokes  equations  in  cartesian  coordinates 


P 
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(2.15) 


where  x  is  the  rate  of  strain  tensor 
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2^u  du  ^  ^  ^  ^  ^ 
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For  the  three-dimensional  equations,  the  following  set  of  nondimensional  variables 

are  introduced 
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which  are  partially  different  from  the  two-dimensional  set  of  the  previous  section. 

The  vorticity  field  is  defined  as  ©  =  -  VxU,  and  with  the  above  nondimensionaliza- 

tion  the  vorticity  components  are 

_ 

*  Re  dz  ^  ' 

(2.18a) 

dw  du 

(2.18b) 

_  du _ 1  dv 

dy  Re  dx  ‘ 

(2.18c) 

By  taking  the  curl  of  the  momentum  equation  (Z15)  and  exploiting  the  foct  that  for 
incompressible  flows  V-u  =  0,  V-©  =  V-Vxu  =  0  and  V^a  =  -  Vx©,  the  vorticity 
transport  equation  for  a  fluid  with  temperature  dependent  viscosity  is 
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doo  .  1 
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where  x  is  the  nondimensional  rate  of  strain  tensor 
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The  nondimensional  operators  Vi  and  V^  are  defined  as: 
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For  the  calculation  of  the  velocity  components  u,v,w 

the  three  partial 

differential  equations 
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are  used,  where  V^  is  the  two-dimensional  Laplace  operator  defined  by 

^"dx^  dx^' 


(2.23) 


These  Poisson  type  equations  are  derived  using  the  vorticity  equations  (2.18)  and  the 
continuity  equation  V  u  -  0.  This  particular  form  with  u  and  w  treated  differently 
from  the  v  component  is  advantageous  for  the  numerical  method  employed  (Fasel  et. 
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aL  1987).  Using  an  implicit  numerical  scheme,  the  Poisson  equations  for  the  stream- 
wise  velocity  component  u  and  the  spanwise  component  w  can  be  solved  directly. 

Hnally  the  thermal  energy  equation  in  nondimensional  form  is 

+  u  -V  e  =  0  ,  (2.24) 

at 

where  6  is  the  nondimensional  temperature,  and  the  viscous  dissipation  has  been 
dropped  as  in  the  two-dimensional  equation,  since  the  Eckert  number  is  very  small. 

2.Z1  Botmdary  Conditions 

The  numerical  simulations  presented  in  this  work  are  based  on  the  assumption 
of  a  two-dimensional  base  flow  subject  to  three-dimensional  disturbances.  These 
three-dimensional  disturbances  are  assumed  to  be  spanwise  periodic.  Then  the  boun¬ 
dary  conditions  for  the  steady  flow  are  the  same  as  the  boundary  conditions  dis¬ 
cussed  in  section  (2.1.2). 

At  the  plate  surface,  A-B  in  figure  2,  the  conditions 
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are  used.  At  the  plate  surface  two  equations  for  the  streamwise  and  the  spanwise 
wall  vortidty  have  to  be  solved.  A  discussion  of  the  choice  of  conditions  for  the  wall 
vortidty  components  can  be  found  in  Easel  et.  al.  (1987).  The  generation  of  three- 
dimensional  waves  in  the  boundary  layer  with  heater  strips  is  discussed  in  the  next 
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chapter. 

At  the  inflow  boundary,  A-D  in  Figure  2,  the  same  conditions  as  for  the  purely 
two-dimensional  flow  are  used.  The  three-dimensional  disturbances  are  introduced  in 
the  flow  field  downstream  of  the  inflow  boundary  ming  heater  strips.  Then  all  dis¬ 
turbance  flow  quantities  are  zero  at  the  inflow  boimdary,  and  for  the  base  flow  quan¬ 
tities  the  similarity  (Blashis)  solution  is  used. 

At  the  top  of  the  integration  domain,  C-D  in  Rgure  2,  exponential  decay  for  all 
disturbance  components  of  the  velocity  is  enforced  similar  to  the  two-dimensional 
Ccise.  Also  the  flow  is  assumed  irrotational,  so  that  all  vorticity  components  vanish. 
The  following  conditions  are  used. 
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w'  ,  (2.26a,b,c) 
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where  the  prime  quantities  denote  dishuiTance  flow.  The  choice  of  g*  will  be  dis¬ 
cussed  later  in  the  context  of  the  numerical  method. 


At  the  outflow  boundary,  B-C  in  Figure  2,  the  conditions  for  the  disturbance 
components  of  the  velocity  «md  the  vorticity  are 
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The  outflow  boundary  conditions  (Z27)  are  chosen  so  that  no  disturbance  com¬ 
ponents  reach  the  downstream  boundary.  A  special  numerical  treatment,  which  will 
be  discussed  in  the  next  chapter,  ensures  that  no  adverse  effects  are  introduced 
•  through  these  conditions  during  the  numerical  integration. 
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3.  Numerical  Method 

For  the  solution  of  the  set  of  partial  differential  equations  given  in  the  previous 
chapter,  two  computer  codes  were  devetoped,  one  for  the  solution  of  the  two- 
dimensional  problem,  and  one  for  the  three-dimensional  problem.  For  the  numerical 
solution  of  the  two-dimensional  problem,  a  fully  implicit  finite-difference  method  of 
fourth-order  accuracy  in  streamwise  and  wall  normal  direction  and  second  order 
accuracy  in  time  was  devetoped.  The  basis  for  this  was  a  numerical  method 
devetoped  by  Fasel  (1976)  for  the  solution  of  the  Navier-Stokes  equations  for  the  flat 
plate  botmdary  layer,  and  the  extension  to  fourth-order  accuracy  by  Bestek  (1980)  for 
plane  channel  flow.  The  numerical  method  is  based  on  the  vorticity  transport  equa¬ 
tion  and  two  Poisson  equations  for  the  velocity  components.  Following  their 
approach  for  the  solution  of  the  Navier-Stokes  equations,  the  energy  equation  is 
incorporated  in  the  niunerical  scheme.  Also  the  additional  viscosity  derivative  terms, 
which  all  vanish  in  the  case  of  isothermal  flow,  have  to  be  treated  in  the  vorticity 
transport  equation  (2.3). 

Hrst  a  computer  code  for  second  order  acctiracy  in  streamwise  and  wall  normal 
directions  was  devetoped  based  on  the  equations  2.3,  2.5  and  2.6,  (Dittrich,  1985).  For 
the  extension  to  fourth-order  accuracy,  special  attention  had  to  be  given  to  the  spatial 
finite-difference  approximations  for  the  boimdary  points  and  for  the  grid  points  adja¬ 
cent  to  the  boundary. 

For  the  three-dimensional  computations.  Krai  (1988)  developed  a  numerical 
method  for  non-isothermal  boundary  layer  flows  based  on  the  method  of  Fasel  et.  al. 
(1987).  This  method  uses  a  spectral  decomposition  method  in  span  wise  direction, 
which  is  based  on  the  set  of  equations  introduced  in  chapter  2.2. 
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3.1  Calculation  of  Two-E)iinensional  Flows 

For  the  discretization  of  the  equations  (2.3),  (2.5)  and  (2.6)  the  rectangular 
integration  domain  is  divided  into  finite  intervals  of  length  Ax,  Ay  (Figure  2),  and  the 
time  coordinate  is  discretized  using  the  interval  At.  The  coordinates  of  the  grid 
points  in  the  streamwise  (x),  wall  normal  (y)  and  the  time  (t)  directions  are  specified 
by  the  integers 

X  =  n  Ax  ,  0<n  , 

y  =m  Ay  ,  O^m^M, 
t  =  p  At  ,  p  =  0,1,2  ■  ■  ■  . 

The  left  boundary  of  the  integration  domain  is  at  the  distance  xq,  measured  from  the 
leading  edge  of  the  plate.  The  distance  Xq  is  chosen  so  that  the  inflow  boundary  is  at 
a  position,  where  the  disturbance  waves,  according  to  linear  stability  theory,  pass  into 
a  region  of  amplification.  A  list  of  the  finite-difference  approximations  used  in  the 
numerical  method  is  given  in  Appendix  B  for  the  energy  equation.  The  finite- 
difference  approximations  for  the  Navier-Stokes  equations  and  the  Poisson  equations 
for  fourth-order  accuracy  in  streamwise  and  wall  normal  direction  can  be  found  in 
Bestek  (1980). 

The  discretized  equations  are  solved  using  a  line  iteration  method  with  overre¬ 
laxation.  At  each  line  the  vorticity  transport  and  the  Poisson  equation  for  the  wall 
normal  velocity  component  are  solved  simultaneously.  Then  the  Poisson  equation  for 
the  streamwise  velocity  component  and  the  energy  equation  is  solved.  The  base  flow 
is  calculated  first  using  a  starting  solution  as  obtained  from  the  similarity  solution  of 
the  unheated  boundary  layer,  with  T^  for  the  heater  temperature  of  the  boundary 
condition  (2.10)  set  to  zero.  The  nonlinear  terms  are  upgraded  at  each  iteration  level 
and  the  iteration  for  the  base  flow  solution  is  repeated  imtil  the  difference  of  two 
consecutive  iterations  for  all  the  flow  variables  reach  the  threshold  e<10“*. 
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After  the  base  flow  has  been  calculated  for  a  certain  heater  strip  configuration,  a 
teznperatiu'e  disturbance  is  introduced  into  the  flow  field  through  the  boundary  con¬ 
dition  at  the  wall  modeling  the  heater  strips.  At  each  time  step,  the  integration  of  the 
discretized  equations  is  similar  to  that  for  obtaining  the  steady  solution.  However, 
now  at  each  time  step  the  iteration  continues  imtil  the  threshold  of  e<10~^  for  the  all 
flow  variables  is  reached. 

For  the  wave  packet  simulations  the  number  of  time  steps  to  be  computed  is 
determined  by  the  speed  of  propagation  of  the  wave  packet  disturbances.  The  max¬ 
imum  number  of  time  steps  is  chosen  so  that  at  in  the  computational  domain  the  last 
X  station  of  interest,  the  disturbances  have  completely  passed  that  position. 

3.2  Calculation  of  Three-CMmensional  Flows 

For  the  three-dimensional  computations  a  spanwise  periodicity  of  the  distur¬ 
bance  components  is  assumed.  Then  the  three-dimensional  equations  can  be  spec¬ 
trally  decomposed  in  spanwise  direction  by  a  Fourier  series 

k^K 

f(pc,y;z,t)=  Z  fkOt,y,f)  ,  (3.1) 

k^K 

where  f  (x,y,z,t  )  represent  all  the  dependent  variables  in  the  physical  domain, 
Fic(x,}/ 4)  are  the  complex  spectral  components  and  y  is  the  spanwise  wavenumber. 

With  the  above  Fourier  approximation  the  vorticity  transport  equation,  the 
energy  equation  and  the  Poisson  equations  for  the  velocity  components  are 
transformed.  Then  all  derivatives  with  respect  to  the  spanwise  direction  are  replaced 
by 

=  iyk  fk(x,y,f)  ,  (3.2a) 

az^  k—K 


(3.2b) 
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The  nonlinear  terms  in  the  vortidty  transport  equations  require  a  special  attention. 
For  a  series  truncation  with  a  large  number  K  each  quadratic  term  in  the  equations 
would  require  a  expansion  proportional  to  terms  to  fully  accotmt  for  the  non¬ 
linear  interaction.  Our  computations  shall  be  restricted  to  the  study  of  the  initial 
three-dimensional  stages  where  only  a  small  niimber  of  nonlinear  interactions  are 
relevant  and  therefore  K  can  be  kept  small.  The  nonlinear  terms  are  calculated  with  a 
purely  spectral  method. 

Substitution  of  the  series  representation  (3.1)  into  equations  (2.19),  (2.22)  and 
(2.24)  yields  a  new  set  of  equations  of  the  form 


3Qifc 

3f 

(3.3a) 

3£2yJ. 

3f 

-Ayk'^§.yk'^Q.yk  > 

(3.3b) 

3Q^ 

3f 

—  ^zk'^^zk'^Q.Af 

(3.3c) 

dx^ 

(3.3d) 

'^2  V 

(3.3e) 

3x2 

(3.3f) 

30^ 

3f 

(3.3g) 

and 

Re  3*2  3y2  Re  ' 

In  the  equations  above  all  variables  are  spectral  coefficients  as  defined  in  (3.1).  The 
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term  ^  is  a  sum  of  all  the  nonlinear  terms  of  the  vorticity  transport  equation,  ^ 
accounts  for  the  diffusion  term  which  is  multiplied  by  the  viscosity  coefficients,  and 
Cji  combines  all  the  additional  viscosity  derivative  terms  in  the  vorticity  equation.  The 
term  Dj  in  the  energy  equation  includes  all  the  convective  terms.  With  this  spectral 
approximation  (3.1),  the  a  required  in  the  boimdary  conditions  for  the  free-stream 
boimdary  conditions  (2.26)  is  a*  =  (c^  +  A  complete  list  of  the  terms 

and  Dj  can  be  foimd  in  Krai  (1988). 

The  three-dimensional  temperature  disturbances  &  are  assumed  to  be  sym¬ 
metric  with  respect  to  the  z=0  axis.  Then  the  disturbance  variables 
symmetric  to  the  z=0  axis  in  the  spanwise  direction.  Since  all  variables  /(x,y,z,f)  in 
(3.1)  are  real,  then  the  U',  V',  O'*  are  purely  real,  while  Q',  ,  £2'y  and  W  are  purely 
imagiruiry.  For  the  complete  solution  of  the  system  of  equations  (3.3)  it  is  therefore 
only  necessary  to  compute  the  spectral  coefficients  from  k=0  to  k=K,  since  for  a  real 
valued  function  f,  the  spectral  coefficients  follow  the  relation  F_jfc  =F*jt,  where  the 
asterisk  denotes  the  complex  conjugate. 

Once  the  equations  are  spectrally  decomposed  in  the  spanwise  direction,  the 
implicit  solution  is  analogous  to  that  for  the  two-dimensional  flow  which  weis  dis¬ 
cussed  previously.  However,  for  the  three-dimensional  calculations,  at  each  grid  line, 
during  the  x-direction  sweep,  K  equations  (3.3)  have  to  be  solved  for  each  spanwise 
coefficient,  before  the  procedure  can  advance  to  the  next  grid  line.  Thus,  for  the 
three-dimensional  calculations  the  solution  has  to  completely  solve  for  a  N  x  M  x 
(K+1)  grid,  before  proceeding  to  the  next  iteration  level. 

For  the  control  of  three-dimensional  disturbances,  the  boundary  conditions 
modeling  the  heater  strips  have  to  incorporate  the  spanwise  variation.  For  the  present 
investigations  the  boundary  conditions  for  the  wall  temperature  (2.25a)  are 
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^wM.k^  =  (  l+^2D^c  ) 


3tOc-Ji) 

X2-X1 


(3.4a) 


^wttU,  k*i  ®  ^3D  Tc  sin 


iiCc-Xi) 

X2-X1 


(3.4b) 


where  the  heater  strip  location  is  in  the  interval  (xi ,  X2],  and  ^420  and  are  the 
amplitudes  of  the  two-  and  three-dimensional  disturbances,  respectively.  For  the 
evaluation  of  a  control  strategy  to  attenuate  three-dimensional  wave  packet  distur¬ 
bances,  a  flow  perturbation  was  triggered  using  only  one  spanwise  mode.  Thus  for 
such  a  control  strategy  to  be  practical  one  would  have  to  assume  that  dominemt 
three-dimensional  modes  are  present  in  the  flow  field.  Experimental  investigations 
have  shown  that  at  initial  three-dimensional  stages  of  transition,  preferred  spanwise 
modes  are  indeed  present.  However,  for  a  more  advanced  flow  control  strategy, 
applicable  to  later  stages  of  the  transition  process,  several  spanwise  modes  would  be 
required. 


The  outflow  condition  at  the  boundary  B-C  in  Figure  2  poses  a  special  problem 
in  the  computation,  since  it  is  actually  not  a  physical  boimdary.  In  order  to  avoid 
that  disturbance  components  reach  the  outflow  boundary  and  become  reflected,  a 
similar  numerical  procedime  was  applied  as  discussed  in  Kloker,  Konzelmann  and 
Fasel  (1991).  At  each  iteration  level  a  tapering  function  is  imposed  for  the  spanwise 
disturbance  vorticity  oOz,  so  that  in  the  last  10%  of  the  integration  domain,  the  span- 
wise  disturbance  vorticity  is  forced  to  zero.  The  tapering  function  used  for  this  pro¬ 
cess,  is 


®'z  k  =  ®'z  k(l  -^Ox^  -6x5)  (35) 

With  this  condition  imposed,  also  all  the  other  disturbance  variables  are  diminished 
within  the  last  10%  of  the  integration  domain.  Several  tests  for  this  boundary  condi- 


_  I 
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tion  have  been  performed.  It  was  confirmed,  that  the  upstream  effect  of  this  tapering 
filiation  remains  confined  to  the  last  10%  of  the  integration  domain,  as  required. 
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4.  Results 


4.1  Computation  of  Two-Dimensional  Flows 

The  first  part  of  the  numerical  results  considers  a  two-dimensional  base  flow 
with  two-dimensional  disturbances.  In  this  case  the  set  of  equations  is  reduced  to  a 
scalar  vorticity  transport  equation  for  the  spanwise  vortidty  . 

For  the  two-dimensional  computations,  the  following  parameters  have  been 
used: 


Re  =  100000 , 
Re^iat  xq)  =  600 , 
At  =  1.122x10-2 , 
Rejjfaf  ofjv)  =  920 , 


L  =  .05  m , 

Ax  =  8.267x10-2 , 

a  =  38 , 

T«  =  20®C  , 

Y/Ay  =  80 , 


xo  =  1.2169 , 

Ay  =  Z1345X10"'  , 
T„  =  30°C  , 

Pr  =  7 
X/Ax  =200  . 


Figure  3  shows  the  location  of  the  integration  domain  in  the  linear  stability 


diagram  frequency  parameter  F  versus  Re^.  The  nondimensional  frequency  parame¬ 


ter  F  is  defined  as 


F=l|=-10^,  (4.1) 

where  P  is  the  dimensional  angular  frequency.  However,  the  frequency  of  the 
Navier-Stokes  simulations  is  nondimensionalized  by 

P  =  (4.2) 

Then,  the  relation  between  these  two  frequencies  is  given  by 

(«) 


where  Re  is  the  Reynolds  number  based  on  a  characteristic  length  L. 
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Based  on  wavenumber  of  the  most  unstable  frequency  of  a  wave  packet  dishir- 
bance  for  the  Res,  range,  covered  in  this  computation,  the  streamwise  stepsize  Ax 
was  calculated  from 


Ax  = 


2ic 


aN, 


(4.4) 


with  NgMw  =  20  for  sufficient  streamwise  resolution  of  the  disturbance  components. 
Therefore,  the  integration  domain  spans  approximately  ten  wavelengths  in  the 
streamwise  direction.  The  height  of  the  integration  domain  Y  was  calculated  from 
the  displacement  thickness  at  the  inflow  boundary  so  that  Y  =  9  5]. 


For  <dl  simulations,  two  heater  strips  were  used  on  the  plate  surface.  The  first 
heater  strip  is  strictly  used  to  create  the  initial  wave  packet  disturbances.  The  second 
heater  strip,  placed  further  downstream  on  the  flat  plate,  is  used  to  control  the  dis¬ 
turbances.  For  the  control  aspects  discussed  in  this  work,  the  creation  of  the  initial 
wave  packet  in  the  flow  field  could  have  also  been  achieved  by  other  methods,  (for 
example  suction-blowing  through  wall  slots).  However,  aspects  of  the  initial  creation 
of  wave  packet  disturbances  with  heaters  could  also  have  an  effect  on  the  control 
implementation. 


4.1.1  Wave  Packet  Disturbances 

In  a  first  attempt  to  create  wave  packet  disturbances,  the  first  heater  strip  was 
activated  by 


=  Tc(f)  sin^ 

X2-X1 


Tc  = 


TP 


X]  ^  X  ^  X2  , 


0  it  TP 
t>TP  , 


(4.5) 


where  TP  is  the  duration  of  the  temperature  forcing  of  the  heater  strip. 
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Rgure  4  shows  the  vorticity  disturbance  signal  at  the  wall  over  the  whole 
integration  domain  at  several  time  steps  after  the  forcing  is  completed.  The  heater 
strip  is  located  between  20  ^  x/dx  ^  30.  A  wave  packet  disturbance  can  be  identified, 
which  travels  in  streamwise  direction.  Fiuthermore,  besides  this  ToUmien  Schlichting 
wave  mechanism,  a  second  effect  is  clearly  identifyable.  A  strong  vorticity  defect  at 
an  early  time  in  the  vicinity  of  the  heater  strip  indicates  a  local  departure  from  the 
initially  isothermal  mean  flow.  The  forcing  signal  (4.5)  introduces  a  local  heat  spot, 
which  is  due  to  the  unsteady  heating.  This  heat  spot  convects  slowly  in  downstream 
direction,  and  diffuses  in  the  wall  normal  direction.  Eventually,  the  spot  disappears, 
and  the  isothermal  condition  reappears  over  the  heater  strip. 

For  the  prospective  control,  the  effect  of  this  unsteady  heating  phenomena  is 
disadvantageous  for  two  reasons.  First,  the  local  non-stationarity  appears  in  the  vorti¬ 
city  signal  as  a  very  low  frequency  band  which  would  have  to  be  filtered  out  by  a 
high  pass  filter  to  identify  the  proper  TS  mechanism.  Second,  the  unsteady  heating 
can  locally  affect  the  stability  characteristics  of  the  mean  flow. 

To  minimize  the  effect  of  such  an  unsteady  component  tHe  heater  strips  are 
operated  over  a  steady  temperature  component.  Then,  a  temporal  disturbance  com¬ 
ponent  is  superimposed  so  that  the  total  mean  heat  input  is  minimized.  The  heating 
procedure  for  the  heater  strips  is  therefore  modified  to  incorporate  a  steady  heater 
temperature  by  using 

e«d,(x,f)  =  (l+Tc(f))sin2(-^^^^^)  ,  (4.6) 

xr-X} 


in  the  interval  (xi  ^  x  ^  X2). 

The  forcing  temperature  (f )  is  taken  to  be 
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24t^ - 45***  +  20t^  if  0<t  ^  TP  12 

Tc  (f )  =  ■  24(2-f  f  -  45(2-f  )*  +  20(2-t  f  if  TP  12  <t  <TP  (4.7) 

0  if  f  >  TP  . 

The  modulus  of  the  forcing  temperature  versus  F  is  shown  on  the  left  hand  side  of 
Bgure  3.  This  particular  choice  ensures  that  the  spectrum  has  a  broad  amplitude  dis¬ 
tribution  over  the  range  of  typical  frequencies  of  the  boundary  layer  and  was  derived 
so  that 

TP 

|T,  (f)df=0  .  (4.8) 

The  effect  of  a  steady  temperature  distribution  on  the  steady  flow  profiles  for 
two  heater  strips  is  shown  in  Figure  5.  A  thermal  wake  behind  the  two  heater  strips 
now  appears  and  reaches  far  downstream.  The  thermal  wake  lies  within  the  boun¬ 
dary  layer  since  for  Pr=7  the  ratio  of  the  thermal  to  the  fluid  boundary  layer  is 
approximately  one  half.  However,  the  magnitude  of  the  temperature  in  the  thermal 
wake  is  small  compared  to  the  temperature  at  the  heater  strip.  Therefore,  the  isother¬ 
mal  contoius  in  Rgure  5  are  chosen  at  a  lower  temperature  range  to  visualize  the 
small  thermal  wake.  The  temperature  gradients  at  stations  further  downstream, 
although  noticeable,  are  small.  The  line  contours  in  Rgure  5  denote  the  constant 
streamwise  velocity.  Over  the  heater  strip  a  slight  indentation  of  the  streamwise  velo¬ 
city  is  apparent,  however,  the  change  of  the  velocity  profiles  downstream  of  the 
heater  strips  is  very  small. 

Results  of  a  computation  with  the  revised  heating  procedure  (4.7)  are  shown  in 
Rgure  6.  The  heater  strip  location  is  between  30  ^  r/Ax  ^  40.  The  wave  packet  distur¬ 
bance  propagates  in  the  downstream  direction,  and  the  local  effect  of  the  unsteady 
heating  is  greatly  reduced.  Only  at  time  f /Af  =  80  can  a  small  non-stationary  heat 
spot  be  identified. 
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A  more  comprehensive  look  at  the  evolution  of  the  wave  packet  disturbance  is 
given  in  Figure  7.  The  disturbance  components  are  plotted  in  perspective  representa¬ 
tion  for  five  time  instances,  f/A/=40,  80, 120,  160,  200,  (the  forcing  ended  at  f/A/=40). 
The  stream  wise  disturbance  velocity  component  at  approximately  y/5  =  0.4  (5  is 
the  local  boundary  layer  thickness),  shows  the  typical  group  of  disturbance  extrema 
in  the  vicinity  of  the  wall,  and  a  smaller,  second  group  of  disturbance  extrema 
further  out,  where  y/5  =  1.2.  The  wall  normal  velocity  component  v'  in  comparison 
has  only  one  group  of  disturbance  extrema.  The  vortidty  disturbance  shows  two 
groups  of  disturbance  extrema,  with  the  larger  one  at  the  wall  surface.  At  the  time 
t/Af  =40,  for  the  u'  velocity  and  for  the  vortidty  signal  co'^  a  spike  close  to  the 
heater  strip  is  visible.  At  the  same  time,  the  disturbance  temperature  O'  shows  two 
temperature  spikes,  which  are  then  strongly  reduced  and  eventually  disappear  over 
the  next  four  time  levels.  Then  the  temperature  wave  packet  emerges  at  a  very  low 
amplitude  level  which  propagates  in  streamwise  direction. 

With  the  revised  temperature  forcing  procedure,  a  slight  deviation  from  a  sta¬ 
tionary  base  flow  is  still  present.  However,  the  TS  wave  mechanism  dominates 
immediately  after  the  end  of  the  forcing. 

From  a  complete  simulation  one  can  analyze  the  wave  packet  disturbances  by 
decomposing  the  signals  using  a  fast  Fourier  transform  (FFT).  The  details  of  using 
such  a  FFT  on  non-periodic  signals  are  explained  in  Appendix  A. 

For  example,  the  wall  vortidty  disturbance  signal  from  the  computed  data  at 
several  streamwise  positions  were  Fourier  analyzed.  Figure  8  shows  the  amplitude- 
and  phase-spectra  of  the  wall  vortidty  when  only  the  first  heater  was  activated.  The 
amplitude-  and  phase-spectra  were  computed  at  several  streamwise  positions  x/Ax  = 
50,  80,  110,  140,  170.  In  the  lower  frequency  range,  for  F  ^  .8,  and  for  frequendes 
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above  Fs  2.0,  each  individtial  frequency  component  becoir\es  damped,  since  the  ^g- 
nitude  of  their  amplitudes  for  increasing  streamwise  positions  are  reduced.  In  the 
intermediate  frequency  region  (.8  ^  F  2.0)  the  amplitudes  of  the  frequency  com¬ 
ponents  of  the  wave  packet  disturbance  increase  The  most  amplified  frequency  com¬ 
ponent  shifts  to  lower  frequencies  from  F  =  2.0  at  xl&x=  50  to  F  =  1.4  at  x/Ax  =  170, 
and  the  band  of  amplified  frequency  components  narrows  considerably.  The  phase- 
spectrum  confirms  the  convective  instability  mechanism  through  the  rapid  change  in 
phases  of  the  frequency  components  for  increasing  streamwise  distance  fr’om  the 
heater  strip.  At  a  wall  normal  distance  of  y/Ay  =  20  ( y/5  =  1.2),  the  amplitude-  and 
phase-spectra  from  analyzing  the  wall  normal  velocity  component  is  shown  in  Figure 
9.  The  amplitude-  and  phase-spectra,  which  are  analyzed  at  a  wall  normal  distance 
outside  the  boimdary  layer  y/Ay  =  20,  agree  with  the  results  of  the  wall  vortidty 
analysis. 

In  Figures  8  and  9,  no  single  point  of  neutral  stability  of  the  frequency  com¬ 
ponents  can  be  identified.  Instead,  several  points  of  neutral  amplification  exist,  which 
can  be  defined  by  the  intersection  of  the  amplitude  spectra  for  two  coirsecutive 
streamwise  positions.  These  neutral  amplification  points  are  shifted  to  lower  fre¬ 
quency  components  for  increasing  downstream  positions.  The  non-parallel  boundary 
layer  continuously  alters  its  stability  characteristics  as  can  be  seen  from  the  stream- 
wise  dependence  of  the  maximum  amplification  rate,  and  from  the  streamwise 
dependence  of  the  neutral  points  of  the  frequency  components  in  the  spectra. 

From  the  phase-spectrum  0  one  can  extract  information  on  the  dispersive  nature 
of  the  wave  packet  evolution.  The  wavenumber  a  of  each  frequency  component  can 
be  foimd  by 


®  Bx  ■ 


(4.9) 
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The  phase  velocities  then  follow  from 


(4.10) 


The  computed  dispersion  relation  a  =  a(p)  of  the  wall  vortidty  signal  at  the 
position  x/Ax  =  140  ,  Reg]  =  835  is  shown  in  Figure  10.  The  range  of  the  curve  is  lim¬ 
ited  by  the  band  of  detectable  frequencies  in  the  wave  packet  spectrum.  For  the 
range  of  frequencies  that  constitute  the  wave  packet  disturbance,  the  dispersion  rela¬ 
tion  is  not  exactly  a  straight  line.  However,  the  departure  from  a  straight  line  is 
smalL 


In  Hgure  11,  a  comparison  is  made  of  the  phase  velocities  and  the  amplification 
rates  obtained  from  the  Navier-Stokes  computations  at  a  position  x/Ax  =140  with 
those  of  the  linear  stability  theory.  The  amplification  rate  a,,  which  denotes  the  ima¬ 
ginary  part  of  the  complex  wavenumber,  is  calculated  from 

a,=^Iny4,  (4.11) 

dx 

where  A  are  the  amplitudes  of  the  frequency  components.  The  phase  velocities  and 
the  amplification  rates  in  Figure  11  were  determined  by  analyzing  the  wall  vorticity 
signal  of  the  Navier-Stokes  equations.  For  the  range  of  frequencies  present  in  the 
wave  packet,  the  phase  velocities  of  the  linear  theory  and  of  the  computations  agree 
very  welL  The  rate  of  amplification  in  Figure  11b  are  identical  for  almost  all  fre¬ 
quency  components,  which  can  be  detected  in  the  wave  packet.  However,  at  the 
maximum  rate  of  amplification  of  the  linear  theory  is  slightly  under  the  value  found 
from  the  Navier-Stokes  computation. 

In  Hgures  12  and  13,  a  comparison  of  the  amplification  rates  of  the  Navier- 
Stokes  computation  at  the  location  of  the  second  maximum  of  the  vorticity  and  of 
the  maximum  of  the  wall  normal  velocity  with  the  linear  theory  is  presented.  The 
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amplification  rates  of  the  wall  normal  velocity  component  exceeds  the  linear  theory 
slightly,  and  also  lie  above  the  prediction,  from  the  wall  vorticity  analysis.  However, 
the  Navier-Stokes  solution  of  the  vorticity  at  the  second  maximum  shows  a  strong 
underprediction  of  the  amplification  rates  compared  to  the  linear  stability  theory. 

For  all  flow  quantities  that  were  cinalyzed  at  three  different  locations  in  the  flow 
field,  the  phase  velocities  are  in  very  good  agreement  with  the  linear  theory  (Figures 
11a,  12a,  13a).  However,  the  rate  of  amplification  can  vary  considerably  depending  on 
the  criteria  used. 

The  amplification  rates  of  the  linear  spatial  theory  have  been  calculated  with  the 
Orr-Sommerfeld  equation  (1.0)  using  the  computed  Navier-Stokes  velocity  profile  as 
the  base  flow.  This  profile  includes  the  small  deviation  from  the  imheated  Blasius 
profile  because  of  the  presence  of  the  thermal  wake.  Although  this  effect  at  first 
seemed  negligible,  it  causes  a  difference  of  almost  15%  in  the  amplification  rates  of 
the  most  amplified  components. 

The  predicament  is  choosing  the  proper  flow  component  at  a  certain  location 
with  which  to  compare  the  amplification  rates  with  linear  stability  theory.  An  addi¬ 
tional  source  of  disparity  between  these  results  is  the  assiimption  of  a  locally  parallel 
base  flow  in  the  linear  theory.  Nonparailel  effects  are  strongest  close  to  the  wall,  and 
the  growth  rates  strongly  depend  on  the  criteria  used  (Fasel  and  Konzelmann,  1990). 
As  for  the  computation  of  wave  packet  disturbances  with  a  broad  band  of  frequen¬ 
cies,  the  analysis  of  the  amplification  rates  show  the  same  trend  as  the  findings  of 
Fasel  and  Konzelmann. 

4.1.2  Passive  Heated  Segments 

The  comparison  of  the  computed  wave  packet  disturbances  with  the  linear  sta¬ 
bility  theory  showed  a  significant  deviation  in  the  amplification  rates  when  the 
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Blasius  pro£Qe  instead  of  the  Naviar-Stokes  profile  was  used  in  the  linear  stability  cal¬ 
culations.  This  effect  can  be  attributed  to  the  thermal  wake  downstream  of  the  heater 
strip. 

To  further  investigate  the  effect  of  passive  heating  segments  on  the  stability  of 
wave  packet  distiubances,  a  heated  segment  was  used  as  a  boundary  condition  as 
shown  in  Figure  14.  The  wall  temperature  distribution  is  shown  on  the  lower  part  of 
Figure  14a,  where  the  heated  segment  starts  at  x/Ar  =  70  and  continues  to  the  end  of 
the  integration  domain.  From  the  Fourier  analysis  of  the  wave  packet  signal  of  the 
wall  vorticity,  several  frequency  components  were  chosen  and  their  normalized 
streamwise  amplification  was  plotted  (Rgure  14a). 

For  the  computation  with  the  heated  segment,  all  the  frequency  components  are 
damped  as  soon  as  they  have  passed  the  temperature  ramp  (Figure  14a).  This  agrees 
favorably  with  the  passive  control  computations  of  Bestek  et.  al  (1987).  Heating  of  the 
boimdary  layer  in  water  stabilizes  the  flow  since  the  viscosity  is  reduced  and  the 
velocity  profiles  become  fuller.  However,  during  the  passage  of  the  temperature 
ramp,  where  d  Q/dx  >  0,  the  frequency  components  from  analyzing  the  wall  vorti- 
dty  shoot  up  to  a  higher  level.  Then  they  all  become  damped  shortly  behind  the  tem¬ 
perature  ramp,  where  d  Q/dx  la*//  =  0. 

In  Figure  14b  a  similar  computation  has  been  repeated,  however  without  the 
heated  segment.  In  this  case,  three  frequency  regions  of  amplification  of  wave  com¬ 
ponents  can  be  identified.  First,  the  lower  frequency  range  consists  of  modes  whose 
amplitudes  initially  decrease  in  streamwise  direction  until  a  minimum  is  reached, 
from  which  point  on  the  modes  are  amplified  (see  e.g.  F=.93).  These  frequencies  pass 
through  the  lower  branch  of  the  neutral  stability  loop.  At  higher  frequencies,  certain 
modes  are  first  amplified  as  they  propagate  downstream,  and  then  become  damped 
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(see  e.g.  F=1.89),  which  means  that  they  passed  through  the  upper  branch  of  the  neu¬ 
tral  loop.  Frequency  components  smaller  than  F=.55  are  damped  over  the  entire 
computational  domain,  but  their  damping  rate  is  progressively  reduced  in  streamwise 
direction. 

To  further  study  the  effect  of  the  temperature  ramp,  a  simulation  with  a  passive 
heater  strip  was  investigated.  Figure  15  shows  a  comparison  of  the  streamwise 
amplification  of  four  frequency  components  from  the  wave  packet  disturbance  with 
and  without  a  passive  heater  strip.  The  heater  strip  is  between  70  <  x/Ax  <  90.  The 
frequency  components  from  the  wall  vorticity  (Figure  15  a)  show  the  selective 
amplification  as  described  above.  Furthermore,  over  the  first  half  of  the  passive 
heater  strip,  the  amplitude  of  each  frequency  component  increases  approximately  by 
10%.  Then,  in  the  second  part,  the  amplitude  of  each  frequency  component  drops 
below  the  value  of  the  unheated  case. 

A  comparison  of  the  same  frequency  components  by  analyzing  the  wall  normal 
velocity  component  away  from  the  wall  (  y/Ay  =  10,  Figure  15b)  does  not  show  the 
local  frequency  "hump"  of  the  wall  vorticity.  However,  both  analyses  show  a  stabiliz¬ 
ing  trend  downstream  of  the  heater  strip  that  continues  until  the  end  of  the  integra¬ 
tion  domain.  It  appeeirs  that  the  additional  heating  from  the  passive  heater  strip  rein¬ 
forces  the  thermal  wake.  The  effect  of  this  reinforcement  reduces  the  rate  of 
amplification  of  the  most  amplified  mode  by  approximately  27%  (Figure  16).  The  dis¬ 
tance  where  the  wall  vorticity  sigrul  was  recorded  for  the  computation  of  the 
amplification  rate  is  two  and  a  half  heater  lengths  downstream  of  the  passive  heater. 
As  seen  in  Kgure  16,  the  dome-shaped  amplification  curve  for  all  frequency  com¬ 
ponents  at  the  position  x/Ax  =  140  is  within  the  one  computed  without  the  passive 
heater.  The  additional  passive  heater  therefore  leads  to  an  overall  less  unstable  boun- 
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dary  layer  due  to  the  thermal  wake. 


Figure  17  shows  the  streamwise  velocity  profQe  of  the  base  flow  and  the  base 
flow  vortidty  at  a  position  x /Ax  =100  close  to  the  passive  heater  strip.  Also  shown  in 
this  figure  are  the  differences  of  the  quantities  between  the  case  with  the  passive 
heater  and  the  case  without  the  passive  heater.  A  fuller  velocity  profile  throughout 
the  boundary  layer  can  be  observed  for  the  heated  case.  However,  the  differences  in 
the  velocity  profiles  are  very  small.  Also,  an  inflection  point  in  the  streamwise  velo¬ 
city  in  the  vidnity  of  the  wall  can  be  identified,  where  doa/dy  =  0.  This  inflection 
point  can  only  be  attributed  to  the  thermal  wake,  which  makes  the  fluid  hotter 
within  a  region  inside  the  boundary  layer.  Since  the  momentum  and  energy  equa- 
tioi\s  are  coupled,  the  viscosity  in  the  thermal  wake  is  reduced,  which  contributes  to 
a  fuller  velocity  profile. 


In  Figure  18,  the  steady  streamwise  pressure  gradient  at  the  wall  surfece  for  two 
passive  heater  strips  is  presented.  The  pressure  gradient  at  y  =  0  is  calculated  from 


dp  ,^d^u  ,dv.fdOy,fdu. 


(4.12) 


In  the  absence  of  the  heaters,  the  pressure  gradient  dp/dx  vanishes,  since  no  external 
pressure  gradient  is  imposed  on  the  flow  field.  In  the  vicinity  of  the  heater  however, 
the  sudden  temperature  rise  leads  to  a  negative  d  0/dy  locally,  so  that  the  wall  pres¬ 
sure  rises.  Note  that  dv/dO  for  a  small  temperature  difference  is  approximately  a 
negative  constant.  In  the  second  half  of  the  heater,  the  sudden  drop  in  temperature 
then  leads  to  a  region  of  Itigher  fluid  temperature  near  the  wall,  and  the  temperature 
gradient  d  0/dy  changes  its  sign. 


For  the  local  analysis  of  the  wall  vortidty  signal,  the  "hump"  in  the  frequency 
components  in  Figure  15  could  be  explained  by  the  local  rise  and  drop  of  the  pres- 
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Svire  gradient  In  the  first  half  a  positive,  adverse  wall  pressure  gradient  leads  to  the 
strong  growth  of  the  frequency  components,  and  in  the  second  half  the  stronger 
negative,  favorable  pressure  gradient  reduces  the  amplitudes  of  the  frequency  com¬ 
ponents  to  a  level,  which  is  lower  compco'ed  to  the  imheated  case.  Away  from  the 
wall  heater,  for  positions  downstream  of  the  heater  strip,  this  effect  is  diminished. 
However,  the  growth  rates  of  the  frequency  components  downstream  of  the  heater 
strip  are  slightly  reduced,  as  previously  explained  in  connection  with  Figure  16,  since 
the  thermal  wake  changes  the  velocity  profiles  of  the  base  flow  through  the  tempera¬ 
ture  dependence  of  the  kinematic  viscosity. 


4.1.3  Transfer  Function  Concept 

Boundary  layer  transition  is  a  very  complex  phenomenon  and  many  aspects  are 
not  yet  understood.  However,  in  the  early  stage  of  the  transition,  where  small  dis¬ 
turbance  waves  become  amplified,  nonlinear  effects  can  be  considered  small. 
Although  the  mechanism  by  which  the  heater  strip  creates  disturbances  in  the  boun¬ 
dary  layer  may  be  nonlinear,  the  amplitudes  of  the  heat  input  were  such  that  the 
amplitudes  of  the  disturbance  response  were  very  small  initially,  so  that  the 
amplification  process  can  be  considered  linear.  Here  a  relationship  between  an  input 
disturbance  (the  heater  temperature),  and  an  output  fluctuation  (the  wall  vorticity)  is 
sought  The  transfer  function  H  (  F  ),  which  is  defined  through  the  input-output  rela¬ 
tionship,  can  be  foimd  from 


H(F)  = 


Ox,(F) 


(4.13) 


where  (^^^(F)  and  ^te  complex  amplitudes  of  the  wall  vorticity  signal  and 

the  temperature  input,  respectively.  Once  the  transfer  function  is  determined,  it  can 
be  used  to  predict  the  linear  response  of  the  boundary  layer  flow  to  any  heat  input 
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fluctuation. 

Rgure  19  shows  a  schematic  representation  of  the  transfer  function  technique 
applied  between  two  stations  Xj  and  Xq.  On  the  left  side  of  the  Figure  19,  the  input 
temperature  signal  from  the  heater  strip  is  shown  together  with  the  modulus  of  its 
transform.  The  input  temperature  signal  is  superimposed  at  the  center  location  of  the 
heater  strip  at  Xq  according  to  (2.8).  The  right  hand  side  of  Figure  19  shows  schemati¬ 
cally  the  vorticity  response  and  the  modulus  of  its  trai\sform  at  the  station  x^.  The 
numerical  approximation  of  this  transfer  function  assumes  a  relatively  flat  input  spec¬ 
trum  over  all  the  relevant  botmdary  layer  modes.  A  discussion  of  the  definition  and 
validity  of  the  transfer  function  is  given  in  Appendix  A. 

Bgure  20  shows  the  modulus  of  the  transfer  function  for  three  different  heater 
strip  widths.  The  three  widths  are  chosen  relative  to  the  wavelength  of  the  most 
amplified  mode  of  the  wave  packet  in  the  region  of  the  heaters.  The  heater  lengths 
N  =  10,  20,  30  correspond  to  approximately  1/2,  1  and  11/2  wavelengths  of  the 
most  amplified  frequency  component.  In  all  three  cases,  the  same  input  temperature 
disturbance  was  used,  and  the  response  signal  was  recorded  at  the  same  position 
x/Ax  =  70.  The  strongest  vorticity  response  was  found  for  N=10.  The  heater  strip 
with  N=30  shows  a  trough  in  the  region  of  F=1.9.  Qose  to  this  frequency  component, 
the  phases  were  undergoing  a  shift  of  n,  so  that  partial  cancellation  of  neighboring 
modes  traveling  with  approximately  the  same  phase  velocities  occurred.  The  appear¬ 
ance  of  such  troughs  in  the  transfer  function  is  disadvantageous  for  the  control  stra¬ 
tegy.  Thus,  a  heater  strip  length  of  N=10  was  chosen  for  the  numerical  investigation 
of  this  control  concept. 
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4.1.4  Boundary  Layer  Control  of  Two-Dimensional  Wave  Packets 

Using  the  numerical  transit  function  discussed  previously  a  relation  is  esta¬ 
blished  between  the  wall  vortidty  response  at  a  certain  position  downstream  of  the 
heater  aitd  the  forcing  input  temperature  at  the  heater.  The  transfer  function  is  based 
on  a  linear  model  Caster  (1984)  extended  the  transfer  function  model  to  allow  for 
quadratic  interaction  of  frequency  components.  For  the  high  amplitude  forcing  in  the 
experiments,  using  suction  and  blowing  through  a  hole  in  flat  plate,  the  quadratic 
interaction  term  accovmted  for  approximately  5%  of  the  total  response.  Therefore,  in 
the  control  attempt  introduced  here,  the  initial  forcing  temperatures  are  chosen  to  be 
small  enough,  so  that  the  amplitudes  of  the  resulting  disturbances  are  below  the 
threshold  for  the  nonlinear  effects  to  be  important.  As  a  result,  once  the  effect  of  a 
pulse  temperature  disturbance  is  known  through  the  transfer  function,  it  follows  that 
the  effect  of  other  disturbance  forms  is  known  as  well 

For  the  control  of  wave  packet  disturbances  in  the  boundary  layer  two  heater 
strips  were  assumed  on  the  plate  surfoce.  The  strategy  for  the  attenuation  of  the  dis¬ 
turbances  in  the  boundary  layer  is  as  follows:  In  a  first  computation  a  wave  packet 
disturbance  is  created  in  the  flow  field  using  only  the  upstream  heater  strip.  From 
this  computation  the  amplitude-  and  phase-spectra  of  the  wave  packet  disturbance 
are  found  from  the  wall  vortidty  sigiuils  at  all  grid  points  downstream  of  the  first 
heater.  This  computation  provides  information  for  a  flow  when  no  control  is  applied. 
From  a  second  computation,  where  only  the  downstream  heater  strip  is  activated  and 
the  first  heater  strip  is  left  passive,  the  amplitude-  and  phase-spectra  of  the  vortidty 
response  are  determined.  This  second  computation  provides  the  necessary  infoima- 
tion  for  numerical  calculations  of  the  transfer  function  as  explained  in  the  previous 
section.  Then,  for  the  attenuation  of  wave  packet  disturbances,  a  third  computation  is 
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perfonned  with  both  heater  strips  activated.  The  first  heater  strip  creates  a  wave 
packet  disturbance  with  the  same  temperature  input  signal  as  in  the  first  computa¬ 
tion.  However,  the  second  heater  strip  acts  now  as  a  controller  whose  temperature 
input  has  to  be  determined  so  that  the  disturbances  are  attenuated  downstream  of 
the  controller.  The  control  temperature  is  foimd  by  the  relation 

4oCf)  =  ^5,,(F)/f?(F),  (4.14) 

where  )  denotes  the  complex  vorticity  amplitudes  of  the  first  computation,  and 
i?(F)  is  the  transf^  ftmction  from  the  second  computation.  In  the  frequency  domain, 
the  phases  of  dx^CF)  of  each  frequency  component  are  shifted  by  tc,  yielding  the 
desired  phases  of  the  control  temperature  input  Then  the  inverse  Fourier  transform 
of  the  complex  temperature  amplitudes  (4.14)  yields  the  temperature  input  for  the 
controller  strip  in  the  time  domain. 

In  order  to  establish  the  numerical  transfer  function,  the  streamwise  distance 
Xj-xq  of  the  sensor  location  from  the  center  of  the  controller  strip  was  varied.  For  a 
sensor  position  too  close  to  the  controller  strip,  the  transfer  function  spectrum  is 
influenced  by  the  localfy  high  vorticity  spikes,  which  are  created  directly  at  the  strip 
surface.  If  the  sensor  is  moved  too  far  downstream,  the  characteristic  distribution  of 
the  amplitude  spectrum  has  narrowed  considerably  and  the  transfer  function  H  has 
lost  information  at  the  lower  and  the  upper  frequency  components  due  to  the  very 
low  amplitude  level  in  the  spectrum.  Rgure  21  shows  the  modulus  of  the  transfer 
function  at  several  positions  downstream  of  the  controUer  strip.  Downstream  from 
the  position  x/Ax  =  88,  the  moduli  of  the  transfer  functions  nearly  coincide  within  the 
frequency  region  .7  ^  f  i  2.0,  where  the  curves  are  almost  hoitzontal.  At  frequencies 
F  2  2.,  the  moduli  branch  out  quickly  due  to  the  strong  damping  of  the  frequency 
components  for  large  F.  After  several  tests,  the  transfer  function  at  the  station 
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x/Ax  =  92  was  picked  for  the  control  strategy.  Figure  22  shows  the  temperature  con¬ 
trol  signal  for  the  disturtiance  with  the  characteristic  time  delay  and  the  typical  wave 
packet  shape. 

For  a  computation  with  two  heater  strips  the  control  strategies  which  were 
introduced  above  have  been  applied.  The  first  heater  created  a  wave  packet  distur¬ 
bance  iising  the  forcing  signal  of  the  relation  (4.6),  while  the  second  heater  strip 
used  the  control  temperature  shown  in  Figure  22.  Figure  23  shows  the  amplitude 
and  phase  spectra  obtained  from  the  wall  vortidty  signal  at  several  streamwise  posi¬ 
tions  for  the  case  where  control  was  applied.  The  first  amplitude  spectrum  was  deter¬ 
mined  at  station  x/Ax  =  50  which  is  before  the  controller  strip  and  may  be  used  as 
reference  to  compare  with  the  downstream  spectra.  The  second  amplitude  spectrum 
at  X /Ax  =  80  is  exactly  at  the  end  of  the  controller  strip  and  shows  a  strong  increase 
in  amplitude.  However,  this  increase  is  due  to  the  locally  strong  vorticity  at  the 
heater  surface  which  creates  the  antiphase  wave  packet.  Further  downstream  at  the 
positions  x/Ax  =  110, 140,  and  170,  the  amplitudes  are  reduced  by  more  than  a  factor 
of  ten  as  compared  with  the  imcontroUed  case.  The  wave  packet  signal  of  the  wall 
normal  velocity  component  (Hgure  24)  at  y/Ay  =  20  (outside  the  boxmdary  layer) 
shows  the  same  attenuation  trend  of  the  wave  packet.  Thus  the  control  of  the  wave 
packet  disturbances  is  not  just  a  wall  effect.  The  most  imstable  disturbances  in  the 
wave  packet  were  attenuated  by  about  97%  of  their  imcontroUed  levels. 

An  alternative  approach  to  the  above  control  strategy  could  be  based  on  an 
additional  sensor  placed  upstream  of  the  controUer  strip.  Then  the  control  system 
exhibits  a  time  delay,  which  could  be  found  by  relating  the  waU  vorticity  at  the 
upstream  sensor  to  the  waU  vorticity  at  the  downstream  sensor.  In  addition  to  the 
above  described  strategy  based  on  (4.13),  a  time  shift  operator  similar  to  (4.13)  but 
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based  on  the  wail  vortidty  at  the  upstream  sensor  as  the  input  quantity,  would  be 
used.  This  would  accoimt  for  the  time  delay  of  the  dishirbances  between  the 
upstream  cmd  the  downstream  sensor.  Thus,  for  a  realistic  control  strategy,  the  dis¬ 
turbances  are  picked  up  by  the  upstream  sensor.  Then  with  this  signal  the  required 
input  for  the  actuator  (controller)  is  calculated  (using  the  technique  just  discussed) 
and  the  actuator  heater  strip  is  triggered  accordingly.  Therefore,  for  this  control  stra¬ 
tegy,  the  downstream  sensor  is  only  needed  once  (before  the  actual  control  of  transi¬ 
tion)  to  establish  the  transfer  function  (between  the  actuator  and  the  downstream 
sensor)  and  the  required  time  delay  (between  the  upstream  sensor  and  the  down¬ 
stream  sensor).  During  the  actual  control  of  wave  packet  disturbances,  the  down¬ 
stream  sensor  may  be  used  to  check  the  efficiency  of  the  control  technique,  that  is  of 
the  wave  cancellation. 

The  crucial  component  of  this  control  strategy  is  the  control  element  based  on 
(4.14)  between  the  wall  temperature  at  the  actuator  and  the  wall  vortidty  at  the 
downstream  sensor  station.  Therefore,  in  the  present  work  emphasis  has  been  placed 
on  exactly  this  part  of  the  control  strategy. 

4.2  Computation  of  Three-Dimensional  Flows 

In  the  later  stages  of  the  transition  process  the  instability  waves  become  three- 
dimensional,  and  in  the  experiments  preferred  spanwise  structures  can  be  identified. 
Therefore,  the  three-dimensional  computations  described  here  assiune  that  the  distur¬ 
bance  flow  is  periodic  in  the  spanwise  direction.  With  a  spectral  decomposition  in  the 
spanwise  direction,  the  equations  (3.3)  have  to  be  solved  for  the  K  terms  of  the  series 
representation  of  (3.1).  For  the  computations  of  three-dimensional  wave  packet  dis¬ 
turbances,  as  discussed  in  the  following  sections,  the  maximum  spanwise 
waveimmber  is  K=2.  The  three-dimensional  wave  packets  are  created  by  fordng  the 
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two-dimensional  mode  (k=0)  and  the  first  three-dimensional  span  wise  mode  Gc=l). 
The  wave  packet  computations  here  can  be  viewed  as  a  model  to  simulate  the  eauly 
stages  of  transition  up  to  the  point  when  the  disturbances  become  three-dimensional, 
and  where  dominant  three-dimensional  modes  exist. 

For  the  three-dimensional  computations,  the  following  parameters  were  used: 

Re  =  100000  ,  Ax  =  5.0265x10"^ ,  T„  =  44°C  ,  K  =  2  , 

L  =  .05  m  ,  Ay  =  2.8461x10"^ ,  y  =  20  ,  Y/Ay  =  81  , 

Xo  =  1.2169  ,  Af  =  3.74x10“^  ,  T«  =  24°C  ,  X/Ax  =  901  , 

Re6,(  flf  xo)  =  600  ,  a  =  25,  Pr  =  63  Re8,(  at  x^)  =  1304  . 

The  height  of  the  integration  domain  Y  was  chosen  so  that  at  the  inflow  boun¬ 
dary  y  =  125i.  This  results  in  a  boxmdary  layer  thickness  at  the  outflow  boundary  of 
approximately  1/2  Y. 

4.2.1  Linear  Spatial  Stability  Theory 

For  the  selection  of  the  spanwise  wavenumber  y  of  the  three-dimensional  com¬ 
putations,  the  linear  stability  theory  was  used  to  give  an  overview  of  the  linear 
amplification  of  pairs  of  oblique  modes,  their  phase  velocities  and  the  spanwise  to 
streamwise  wavenumber  ratios  for  several  frequeiKies. 

For  linear  spatial  stability  theory  the  wavenumber  in  equation  (1.0)  is  complex 
a  =  a-i-a,-  and  the  circular  frequency  p  is  real  Equation  (1.0)  is  based  on  the  the 
assumption  of  a  two-dimensional  disturbance  wave.  The  three-dimensional  deriva¬ 
tion  of  (1.0)  follows  from  assuming  three-dimensioiral  waves  of  the  form 
^(y)  exp(a3X -H'ftz  -  Paf),  and  assuming  a  parallel,  two-dimensional  base  flow 
Uf^  =  0,  0).  With  the  aid  of  the  Squire's  transformation  (Drazin  and  Reid, 

1981),  a  three-dimensional,  oblique  disturbance  wave  can  then  be  found  from  equa- 


tion  (1.0)  through 


Oj  =  (ot^3  +  7*3)^'^ ,  ot2  Re2  =  03  Re^  ,  ,  (4.15a,b,c) 

03  02 

where  the  subscripts  2  and  3  refer  to  the  two-dimensional  and  the  three-dimensional 
modes,  respectively,  and  where  Re^  is  the  equivalent  three-dimensional  Reynolds 
number.  The  phase  velocity  c  of  equation  (1.0)  must  then  be  replaced  by 
c  =  which  is  the  three-dimensional  component  in  the  stream  wise  direction. 

Since  we  study  pairs  of  oblique  waves,  which  propagate  with  a  wave  angle 
tan~^(Y/  a)  to  the  streamwise  direction  of  equal  and  opposite  sign,  the  amplification 
rate  in  spanwise  direction  Yi  is  assumed  to  be  zero.  Then,  by  specifying  the  spanwise 
wavenumber  y ,  equation  (1.0)  can  be  used  to  calculate  the  eigenvalues  of  the  pairs  of 
oblique  modes. 

Rgure  25  shows  the  rate  of  amplification,  the  phase  velocity  and  the  ratio  of  the 
spanwise  to  streamwise  wavenumber,  which  are  obtained  from  the  linear  stability 
theory.  The  eigenvalues  are  shown  for  two  streamwise  positions  Rej^  =  800,  1200  and 
for  pairs  of  oblique  modes  with  y  =  10,  20,  30,  40,  and  the  two-dimensional  mode 
with  Y  =  0. 

The  rate  of  amplification  for  the  oblique  modes  decreases  with  increasing 
wavenumber  for  both  Reynolds  numbers.  The  maximum  amplification  shifts  to 
lower  fi'equendes  for  increasing  Reynolds  numbers.  For  both  Reynolds  numbers  the 
highest  rate  of  amplification  is  fbtmd  to  belong  to  a  frequency  component  of  the 
two-dimensional  mode  (y  =  0) .  However,  at  Regj  =  1200  for  frequencies  below  F  =  .5, 
the  oblique  mode  with  y  =  10  shows  a  locally  higher  amplification  rate  than  the  two- 
dimensional  mode.  Although  the  maximum  amplification  of  a  frequency  is  connected 
with  the  two-dimensional  mode  (y  =  0),  at  some  other  frequency  a  three-dimensional 
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wave  can  still  be  more  amplified  than  a  two-dimensional  one.  Although  generally 
more  stable  than  two-dimensional  waves,  the  three-dimensional  waves  with  a  small 
wavenumber  y  can  have  similar  growth  rates  as  the  two-dimensional  waves. 

The  streamwise  phase  velocities  are  calculated  from  the  relation  (4.10)  for  each 
frequency  and  oblique  wavemxmber  y.  For  increasing  y,  the  phase  velocities  in  Figure 
25b  become  higher.  Towards  the  lower  frequencies,  the  phase  velocities  decrease  for 
the  two-dimensional  mode.  However,  for  the  higher  spanwise  waveniunbers  y  S  20, 
the  phase  velocities,  after  first  decreasing,  start  to  increase. 

Also,  Figure  25c  shows  that  the  ratio  of  spanwise  to  streamwise  wavenumber 
y  /  a  does  not  change  from  the  position  at  Rejj  =  800  to  Reg^  =  1200.  At  the  first  Rey¬ 
nolds  number,  the  wavenumber  y=20  at  the  frequency  F=1.2  has  a  ratio  of 
y  /  a  =  .67.  At  this  frequency  the  two-dimensional  mode  has  a  maximum 
amplification  (Figure  25a).  At  the  Reynolds  manber  Reg,  =  1200,  where  the  most 
amplified  modes  have  shifted  to  the  lower  frequencies,  the  ratio  changes  to 
y  /  a  =  1.26  at  the  frequency  F  =  0.6, 

For  the  forcing  of  the  pair  of  two  oblique  modes  with  a  wave  packet  distur¬ 
bance,  the  wavenumber  y  =  20  was  selected  for  the  computations.  A  similar  ratio  of 
spanwise  to  streamwise  wavenumber  was  also  reported  in  the  transition  experiments 
of  Saric  and  Thomas  (1983)  for  single  frequency  components  at  amplitude  levels  of 
.3%  to  .4%  based  on  free  stream  velocity.  With  our  wave  packet  simulations  based  on 
initially  low  amplitude  forcing  of  a  wide  band  of  frequency  components,  we  intend 
to  further  investigate  the  role  of  TS  modes  in  the  transition  process. 

4.Z2  Three-Dimensional  Wave  Packet  Disturbances 

Hgure  26  shows  the  instantaneous  wall  temperature  at  f/Af  =  70  for  the  two 
heater  strips.  The  first  hrater  strip  was  activated  with  a  physical  amplitude  of  the 


forcing  temperature  for  the  two-dimensional  and  the  three-dimensional  mode  of 
A2J)  =  A^d  =  20°C.  The  forcing  period  is  TP  =  120  and  the  second  heater  strip  is  pas¬ 
sive  The  temperature  signal  was  recomposed  by  summing  the  series  (3.1),  and  the 
spanwise  extend  of  the  integration  domain  was  chosen  to  comprise  two  spanwise 
wavelengths.  The  strong  three-dimensional  forcing  of  the  wall  temperature  is  visible 
in  the  two  temperature  peaks  which  are  one  spanwise  wavelength  apart. 

The  disturbance  signal  of  the  spanwise  wall  vortidty  at  several  instances  in  time 
is  shown  in  Figure  27  for  the  two-dimensional  co,  and  the  first  three-dimensional 
spanwise  mode  (O,  At  the  first  timestep  tIM  =  180  (the  end  of  the  forcing  is  at 
tIAt  =  120),  the  trailing  side  of  the  wave  packet  illustrates  the  non-stationarity  of  the 
base  flow  due  to  the  unsteady  heating,  which  is  an  effect  disoissed  earlier  in  conjunc¬ 
tion  with  the  computations  for  purely  two-dimensional  flow.  The  three-dimensional 
wave  packet  propagates  in  the  downstream  direction  and  spreads  due  to  dispersion. 
The  wall  normal  velocity  components  Vj^o  and  v^^i  at  a  wall  normal  location 
y/Ay  =  40  outside  the  boimdary  layer  (Bgure  28)  also  confirm  a  traveling  wave 
packet.  At  the  time  level  t/Af  =  180,  the  imsteady  heating  effect  is  not  visible  since 
the  thermal  wake  remains  confined  inside  the  boundary  layer. 

The  spectralty  decomposed  streamwise  velocity  components  Uk=o  and  Ujt=i  are 
displayed  in  Figiue  29  in  the  vertical  (x-y)  plane  at  z=0.  Shortly  after  the  end  of  the 
wave  packet  excitation  at  timestep  130  (the  duration  of  the  excitation  is  until  timestep 
120),  the  streamwise  disturbance  velocity  contoius  indicate  a  wave  packet  for  both 
the  two-dimensional  (k=0),  and  the  three-dimensional  (k=l)  component.  In  this  early 
stage  the  spatial  selection  of  the  unstable  fi^uendes  is  not  yet  fully  completed. 

At  the  later  timestep  1200  (Figure  30),  the  wave  packet  has  spread  considerably. 
The  three-dimensional  component  has  spread  over  less  cydes  than  the  two- 
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dimensional  one.  This  indicates  th  jt  the  difference  of  phase  velocities  between  indivi¬ 
dual  frequency  components  for  the  three-dimensional  mode  is  less  than  the  difference 
of  phase  velocities  for  the  two-dimensional  mode.  Also,  a  third  group  of  disturbance 
extrema  in  the  velocity  contours  of  the  three^iimensional  mode  can  be  identified  in 
the  immediate  vidiuty  of  the  wall 

For  the  spanwise  vorticity  (O*  contours  in  Figure  31,  the  two-  and  the  three- 
dimensional  modes  appear  more  distinct  than  for  the  streamwise  velocity  at  the  same 
time  level  (f /Af=1200).  The  three-dimensional  spanwise  vorticity  mode  has  a  group  of 
disturbance  extrema  at  the  wall  surface.  Away  from  the  wall,  a  second  group  of  dis¬ 
turbance  extrema  of  the  vorticity  along  a  45®  angle  to  the  x-axis  can  be  distinguished. 
The  magnitude  of  the  disturbance  extrema  away  from  the  wall  is  larger  than  that  of 
the  group  of  extrema  on  the  plate  surface.  In  contrast,  the  two-dimensional  spanwise 
vorticity  mode  has  a  group  of  disturbance  extrema  at  the  wall  surface,  whose  magni¬ 
tude  is  larger  than  a  second  group  of  extrema  which  is  located  at  y/Ay=15. 

The  amplitude-  and  phase-spectra  of  the  spanwise  wall  vorticity  (o^  for  the  two- 
dimensional  and  the  three-dimensional  mode  are  given  in  Hgure  32.  The  spectra  are 
recorded  at  several  streamwise  positions  x/Ax  =  90,  100,  110,  120,  130  downstream  of 
the  first  strip.  At  the  station  closest  to  the  activator  strip,  the  spectra  for  both  the 
k=0  and  k=l  component  are  still  very  broad,  and  are  similar  to  the  temperature  input 
spectrum.  However,  further  downstream  the  shape  of  the  spectra  illustrates  the  selec¬ 
tion  mechanism  of  the  boundary  layer.  The  very  high  and  very  low  frequency  com- 
ponenfa  are  strongly  damped  and  the  characteristic  dome-shaped  spectrum  appears, 
where  the  intermediate  frequency  components  are  most  amplified. 

The  amplitude  spectrum  of  the  k=l  component  shows  a  sharp  minimum  at  a 
frequency  of  approximately  F  =  .75.  The  phase-spectrum  in  the  neighborhood  of  this 
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frequency  undergoes  a  sudden  phase  change  of  so  that  modes  in  the  vicinity  of 
this  trough  cancel  each  other.  The  region  of  aunplification  for  the  k=l  spectrum  is 
shifted  to  lower  frequency  components  compared  to  the  k=0  spectrum.  Also,  the 
amplification  rates  are  reduced. 

The  amplitude  spectra  of  the  wall  normal  velocity  component  at  the  same 
streamvdse  position  and  y/Ay  =  40  (Figure  33)  exhibit  a  similar  behavior  as  for  the 
spanwise  vorticity.  Furthermore,  the  three-dimensional  component  does  not  indicate 
the  sharp  trough  in  the  vicinity  of  F  =  .75  .  From  the  position  x/Ax  =110  on,  the 
amplitude  spectra  of  Figures  32  and  33  are  of  the  characteristic  dome-shape  of  wave 
packets  in  the  boundary  layer.  This  position  is  a  distance  of  two  heater  lengths  away 
from  the  center  of  the  heater  strip  in  downstream  direction. 

The  streamwise  development  of  the  amplitude  spectra  over  a  larger  distance  is 
illustrated  in  Figures  34  and  35  for  the  spanwise  wall  vorticity  and  the  wall  normal 
velocity  component.  The  positions  are  at  x/Ax  =  180,  260,  340,  420,  500.  For  the 
spamwise  vorticity,  the  two-dimensional  mode  amplifies  and  selects  its  frequency 
components  according  to  linear  theory.  The  three-dimensional  mode,  however,  shifts 
to  lower  frequency  components  without  increasing  amplification.  Instead,  at  the 
lower  frequency  range  f  ^  .5,  frequency  components  with  an  initially  very  low 
amplitude  become  amplified.  The  three-dimensional  spanv^dse  vorticity  has  its  first 
group  of  disturbance  extrema  at  y/Ay  =  10,  close  to  the  wall  At  the  wall  surface 
itself,  only  a  smaller  group  of  extrema  is  present  (Figure  31).  By  analyzing  the  span- 
wise  vorticity  signal  at  the  position  y/Ay=10  (Fig.  34c)  with  maximum  disturbance 
amplitude,  the  three-dimensional  mode  shows  an  increasing  rate  of  amplification  of 
the  lower  frequency  components.  As  already  pointed  out  by  the  strictly  two- 
dimensional  computations,  the  rate  of  amplification  can  be  determined  at  any  wall 
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nonnal  position  for  any  flow  quantities.  Although  there  is  no  unique  amplification 
rate  using  the  direct  numerical  simulation  to  compare  with  the  linear  parallel  theory, 
comparisons  are  usually  made  from  disturbance  quantities  at  their  maximum  ampli¬ 
tude  location.  From  our  computational  results  we  find  the  best  agreement  with  the 
parallel  theory  also  for  the  location  with  the  maximum  disturbance  amplitude. 

The  amplitude  spectra  for  the  wall  normal  velocity  component  in  Figure  35  is 
calculated  at  y/Ay  =  40  outside  the  boundary  layer  at  the  disturbance  maxima.  The 
amplitude  spectra  indicate  the  same  stability  characteristics  as  for  the  spanwise  vorti- 
dty.  The  three-dimensional  spanwise  mode  exhibits  the  shift  to  lower  frequency  com¬ 
ponents  together  with  an  increase  in  amplification. 

The  amplitude  spectra  of  all  the  analyzed  flow  components  show  a  strong 
growth  of  the  lower  frequency  components  for  the  three-dimensional  mode.  The 
growth  rate  of  these  lower  frequency  components  is  much  higher  than  the  linear 
growth  rates.  In  contrast,  the  two-dimensional  amplitude  spectrum,  as  for  all  the 
analyzed  flow  components,  do  not  indicate  the  growth  of  these  lower  frequency  har¬ 
monics. 

In  Figure  36  a  comparison  of  the  experimental  amplitude  spectra  with  the  com¬ 
puted  spectra  of  the  linear  spatial  stability  theory  from  Caster  and  Grant  (1975b)  is 
presented.  The  spectra  to  the  left  of  Figure  36  are  the  experimentally  determined  dis¬ 
tributions  of  the  two-dimensional  mode  and  a  three-dimensional  mode  for  the 
streamwise  velocity  component,  which  were  measured  at  a  wall  distance  of  y/5  =  1.1 
outside  the  boundary  layer.  To  the  right  of  Figure  36  the  computed  spectra  of  the 
linear  stability  theory  are  shown. 

From  these  figures  one  can  clearly  identify  a  broadening  of  the  amplitude  sp>ec- 
tra  for  three  three-dimensional  modes  towards  the  lower  frequency  components. 
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which  is  not  present  in  the  linear  spectra.  From  the  first  appearance  of  the  low  fre¬ 
quency  "humps"  at  an  early  station,  the  low  frequency  components  eventually  take 
over  the  whole  spectra,  so  that  at  the  last  station  the  most  amplified  frequency  com¬ 
ponents  are  at  a  lower  frequency  range  than  the  linear  theory  could  predict. 

Also,  and  even  more  surprising,  the  two-dimensional  mode  does  not  exhibit  the 
low  frequency  interaction  mechanism  of  the  three-dimensional  mode.  The  two- 
dimensional  spectrum  follows  closely  the  linear  prediction  up  to  the  last  station. 

Our  niunerical  computations  qualitatively  indicate  the  same  interaction  mechan¬ 
ism  for  the  three-dimensional  mode  as  in  the  experimental  investigation  of  Caster 
and  Grant.  However,  a  comparison  of  our  numerical  simrilation  with  the  experi¬ 
ments  can  only  be  qualitative  for  several  reasons.  First,  the  experimental  study 
covered  a  much  larger  streamwise  region.  Second,  the  wave  packet  disturbance  in  the 
experiment  was  created  through  a  single  hole  in  the  plate  surface,  such  that  an  iso¬ 
lated  three-dimensioiral  wave  packet  was  created  rather  than  a  two-dimensional  wave 
packet  with  a  spanwise  variation. 

In  our  computations,  during  the  forcing  of  the  wave  packet  with  the  two- 
dimensioiuil  and  one  three-dimensional  mode,  a  second  spanwise  mode  with  the 
wavenumber  y  =  40  is  created  also.  This  second  spanwise  mode  results  from  a  non¬ 
linear  interaction  of  the  first  spanwise  mode  with  itself  during  the  temperature  forc¬ 
ing  over  the  heater  strip.  However,  the  amplitudes  of  all  frequency  components  of 
this  second  mode  wave  packet  were  small  compared  to  the  two-dimensional  and  first 
three-dimensional  modes.  All  frequency  components  of  this  second  mode  are 
damped  within  a  short  distance  from  the  heater  strip.  For  the  wavenumber  7,  =  40  of 
the  second  spanwise  mode,  the  linear  spatial  theory  predicts  damping.  The  rate  of 
damping  is  further  increased  for  larger  Re  numbers. 
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The  amplification  of  the  lower  frequency  components  for  the  three-dimensional 
spanwise  mode  indicates  the  presence  of  a  subharmonic  amplification  mechanism. 
The  amplification  of  subharmonic  frequencies  was  found  experimentally  to  exist  for 
periodic  disturbances  in  the  early  three-dimensional  stages  of  the  breakdown  process 
(Saric  and  Thomas  (1983),  and  Kachanov  and  Levcherdto  (1984)).  Corke  (1989) 
confirmed  in  an  experimental  investigation  for  single  frequency  disturbances  at  very 
low  initial  amplitudes,  that  the  spectral  broadeiung  of  the  power  spectra  first 
appeared  in  the  subharmoiuc  region.  There  the  breakdown  process  was  absent  of 
higher  order  instability  mechaiusms,  such  as  proposed  in  the  fundamental  break¬ 
down  theories. 

In  Rgure  37,  a  comparison  of  the  phase  velocities  is  made  between  the  two- 
dimensional  mode  and  the  three-dimensional  spanwise  mode  at  x/Ax=480 
(Resj  =  1021)  for  the  wall  normal  velocity  at  y/Ay  =  40.  Above  the  frequency  F  =  .7 
the  three-dimensional  phase  velocities  agree  with  the  linear  theory.  However,  below 
F  =  .5  the  phase  velocities  are  approximately  .36  Lf„.  In  the  low  frequency  range, 
where  the  strong  growth  of  the  amplitudes  of  the  three-dimensional  modes  appeared, 
the  phase  velocities  deviate  from  the  value  of  the  linear  stability  theory  and  level  to 
approximately  the  constant  value.  The  phase  velocities  of  the  two-dimensional  mode 
agree  closely  with  the  linear  theory. 

To  summarize  the  generation  process  of  three-dimensional  wave  packet  distur¬ 
bances  using  heater  strips,  all  the  flow  variables  have  been  recomposed  for  two  span- 
wise  wavelengths  and  are  displayed  in  the  horizontal  (x-z)  plane  in  Figures  38 
through  43.  The  grey-shaded  contours  were  taken  at  the  time  level  t/At  =  70,  the 
period  of  excitation  being  TP  =  120.  Superimposed  are  contours  of  the  wall  tempera¬ 
ture  disturbance,  where  the  zero  contour  frames  the  physical  dimension  of  the  heater 
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strip. 

The  streamwise  and  wall  nonnal  velocity  components  at  their  disturbance 
extrema  show  the  wave  packet  disturbance  aligned  with  the  temperature  forcing  at 
the  same  spanwise  position  (Figures  38  and  39).  E>uring  the  excitation,  there  is  an 
upstream  effect  in  the  velocity  components  in  addition  to  the  propagation  in  down¬ 
stream  direction.  In  comparison,  the  spanwise  wall  vortidty  (Rgure  43)  remains 
more  confined  to  the  temperature  disturbance  during  the  excitation,  although  the 
wave  packet  has  already  moved  slightly  in  the  downstream  direction.  The  spanwise 
vortidty  disturbance  extrema,  as  for  the  streamwise  and  wall  normal  velodty  are 
aligned  with  the  temperature  disturbaiKre  extrema.  The  strong  coherence  between  the 
spanwise  vortidty  component  and  the  wall  temperature  confirms  their  choice  for  the 
estimation  of  the  transfer  function. 

A  different  picture  is  displayed  for  the  purely  three-dimensional  variables 
w',  co'y  (Figures  40  to  42).  There,  the  disturbance  extrema  are  centered  between 
the  wall  temperature  extrema,  and  these  disturbances  vanish  where  the  local  span- 
wise  temperature  gradient  is  zero.  EXuing  the  forcing  of  the  wave  packets,  the 
streamwise  wall  vortidty  cOx  reaches  far  upstream.  Instantaneously,  a  three- 
dimensional  vortidty  field  with  a  large  G)*  component  is  created  in  the  boundary 
layer. 

The  relation  of  the  streamwise  vortidty  co*  and  the  spanwise  vortidty  at  a 
later  time  t/M  =  1200  is  shown  in  Bgure  44.  The  grey-shaded  contours  indicate  the 
streamwise  vortidty  at  the  wall  surface.  The  overlaid  contours  denote  the  spanwise 
wall  vortidty.  In  regions  were  the  spanwise  vortidty  gradients  are  small,  the  stream- 
wise  wall  vortidty  disappears.  On  the  other  hand,  were  the  spanwise  vortidty  gra¬ 
dients  are  large,  the  streamwise  vortidty  has  its  disturbance  maximum.  The  stream- 
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wise  voitidty  disturbance  is  concentrated  along  the  streamwise  direction  between  the 
the  peak  and  valley  positions  as  marked  by  the  spanwise  wall  vortidty. 


4.Z3  Boundary  Layer  Control  of  Three-I>imensional  Disturbances  Waves 

For  the  control  of  three-dimensional  wave  packet  disturbances  the  same  transfer 
function  strategy  as  for  the  two-dimensional  case  (see  4.13)  is  tised.  Since  the  three- 
dimensional  equations  are  already  decomposed  in  the  spanwise  direction,  an  addi¬ 
tional  transfer  function  relating  the  three-dimensional  spanwise  vorticity  to  the  three- 
dimensional  temperature  forcing  can  be  introduced  analogously  to  (4.13); 


(4.16) 


Once  the  effect  of  a  three-dimensional  temperature  input  is  fotmd  using  (4.16),  then  a 
controller  signal  can  be  created  for  the  attenuation  or  the  reinforcement  of  the  distur¬ 
bance  components.  The  procedure  to  find  the  control  temperature  for  disturbance 
attenuation  is  similar  to  the  procedure  which  was  discussed  in  section  (4.1.4.)  in  con¬ 
nection  with  the  purely  two-dimensional  dishirbance  computations.  However,  now 
each  spanwise  component  has  to  be  treated  individually. 


The  moduli  of  the  transfer  function  for  the  two-dimensional  and  the  three- 
dimensional  mode  at  several  streamwise  positions  are  given  in  Figure  45.  The  two- 
dimensional  mode  exhibits  the  same  behavior  as  in  the  purely  two-dimensional  com¬ 
putations.  The  moduli  have  a  constant  amplitude  ratio  within  the  band  of  linearly 
amplified  frequencies.  For  frequency  components  above  F  =  1.8,  the  moduli  of  the 
transfer  function  change  rapidly  with  respect  to  the  streamwise  direction. 

The  moduli  of  the  three-dimensional  transfer  function,  however,  show  a  pro¬ 
nounced  trough  in  the  vicinity  of  the  frequency  F=.75  .  This  effect  could  already  be 
expected  from  the  amplitude  spectra  of  Figure  32.  The  rapid  phase  change  near  the 
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frequency  component  F=  .7  leads  to  partial  cancellation  in  the  vicinity  of  this  fre¬ 
quency,  since  the  phase  velocities  of  the  neighboring  modes  of  the  wall  vortidty 
response  are  very  close  to  each  oAer. 

4^3.1  Control  of  the  Two-Dimensional  and  Three-Dimensional  Span  wise  Modes 

For  the  attenuation  of  the  two-  and  three-dimensional  wave  packet  disturbances 
a  transfer  function  was  used,  which  was  selected  at  the  streamwise  position 
x/Ai  =230.  The  controller  strip  was  placed  between  the  streamwise  position 
130  S  x/A  S  145.  The  x-station  for  the  transfer  function  is  approximately  two  distur¬ 
bance  wavelengths  (of  the  most  amplified  frequency  in  the  wave  packet)  downstream 
from  the  center  of  the  controller  strip. 

Hgures  46  to  49  show  the  amplitude  spectra  for  the  control  of  the  two-  and 
three-dimensirnal  disturbance  components  of  the  spanwise  wall  vorticity  o)  ^  and  the 
wall  normal  velocity  component  v'  at  several  streamwise  positions.  The  amplitude 
spectra  at  the  position  x/Ax  =  120  is  just  upstream  of  the  controller  strip  and  may 
serve  as  a  reference  amplitude. 

For  the  two-dimensional  mode,  a  90%  reduction  of  the  most  amplified  fre¬ 
quency  component  was  achieved.  The  three-dimensional  spanwise  mode  control 
showed  an  approximately  80%  reduction  for  the  spanwise  vorticity.  The  control 
based  on  the  transfer  function  for  the  three-dimensional  spanwise  vorticity  is  not 
quite  as  successful  as  the  control  for  the  two-dimensional,  spanwise  mode.  The 
three-dimensional  vorticity  has  its  primary  disturbance  maximum  away  from  the 
walk  A  more  successful  estimation  of  the  transfer  function  of  the  three-dimeirsional 
mode  might  be  found  away  from  the  wall  or  using  the  streamwise  vorticity. 
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4.Z3.2  Control  of  the  Two-Dimensional  Mode 

Hnalty,  control  of  only  the  two-dimensional  mode  has  been  attempted.  By 
suppressing  the  two-dimensional  mode  in  the  wave  packet  it  was  hoped  to  gain 
insight  on  the  devebpment  of  the  subharmonic  frequency  modes  which  are  present 
in  the  three-dimensional  disturt>ance  development  of  the  uncontrolled  flow. 

Rgure  50  shows  the  spanwise  vorticity  od'z  in  the  horizontal  plane  at  the  plate 
surface  for  the  attenuation  of  the  two-dimensional  mode  only  (grey-shaded  contours) 
at  the  time  1200.  The  spanwise  vorticity  of  the  imcontroUed  flow  is  overlaid  as  line 
contours.  The  suppression  of  the  two-dimensional  mode  advances  the  spanwise  vorti¬ 
city  at  the  valley  location.  At  the  peak  locations  the  maxima  still  coindde  However, 
the  effect  of  trying  to  cancel  only  the  two-dimensional  mode  on  the  lower  frequency 
interactions  can  be  interpreted  better  from  the  amplitude  spectra. 

In  Figure  51,  the  amplitude-  and  phase-spectra  of  the  three-dimensional  mode  of 
the  vorticity  at  y/Ay=10  and  the  wall  normal  velocity  at  y/Ay=40  are  shown  for 
several  streamwise  positions.  Both  flow  components  are  analyzed  at  their  amplitude 
maxima  in  wall  normal  direction.  The  two-dimensional  mode  of  both  flow  quantities 
was  attenuated  similarly  to  the  "two-mode"  control  case.  The  attenuation  of  only  the 
two-dimensional  mode  leaves  the  linear  amplification  mechanism  of  the  three- 
dimensional  mode  intact.  However,  the  strong  amplification  of  the  low  frequency 
components  of  the  three-dimensional  mode,  which  were  present  in  the  uncontrolled 
case,  are  not  identifiable  now  in  the  amplitude  spectrum.  The  attenuation  of  only  the 
two-dimensional,  spanwise  mode  of  the  wave  packet  disturbance  effectively  delayed 
the  onset  of  the  nonlinear,  subharmonic  type  interactions. 
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5.  Conclusions 

A  numerical  investigation  of  the  control  of  the  transition  from  laminar  to  tur¬ 
bulent  flow  in  a  flat  plate  boimdary  layer  flow  was  imdertaken.  The  numerical 
method  is  based  on  the  complete  Navier-Stokes  equations  and  allows  for  the  spatial 
amplification  of  disturbance  waves  in  a  non-parallel  boundary  layer.  For  the  numeri¬ 
cal  method  finite-differences  are  used  in  the  streamwise  and  wall  normal  directions 
for  the  two-  and  three-dimensional  computations,  and  a  spectral  Fourier  approxima¬ 
tion  is  used  in  the  spanwise  direction  for  the  three-dimensional  computations. 
Emphasis  is  placed  on  the  simulation  of  two-  and  three-dimensional  wave  packet  dis¬ 
turbances  since  they  more  realistically  model  the  natural  transition  process  than  sin¬ 
gle  frequency  disturbances.  In  the  numerical  simulation  flush  mounted  heater  strips 
on  the  plate  surface  were  modeled  to  generate  disturbance  waves. 

We  investigated  the  mechanism  by  which  Tollmien  Schlichting  (TS)  waves  are 
created  in  the  boimdary  layer  through  a  local  viscosity  variation.  The  time  dependent 
heat  transfer  during  the  thermal  forcing  with  the  heaters  locally  alters  the  viscosity, 
which  then  imposes  a  temporal  fluctuation  on  the  streamwise  velocity  profile.  Strong 
vorticity  gradients  over  the  heater  strip  appear  instantaneously,  thereby  generating 
the  TS  waves. 

During  the  temperature  forcing,  in  addition  to  the  TS  mechanism,  a  net  amount 
of  heat  is  introduced  in  the  flow  field.  At  the  end  of  the  forcing,  a  region  of  locally 
high  temperature  (heat  spot)  appears  in  the  vicinity  of  the  heater  strip,  which  is  con¬ 
verted  in  the  streamwise  direction  at  a  speed  much  slower  than  the  group  velocity  of 
the  TS  wave  packets.  The  heat  spot  rapidly  diffuses  in  wall  normal  direction  and 
almost  vanishes  at  a  downstream  distance  of  approximately  two  heater  strip  lengths. 
The  effect  of  this  non-stationarity  appears  in  the  Fourier  analyzed  signals  as  a  low 
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frequency  component,  which  adversdy  influences  our  control  strategy. 

A  steady  temperature  component  over  the  heater  strip  reduces  the  effect  of  this 
non-stationarity  in  the  sensor  signal  However,  this  introduces  a  thermal  wake  in  the 
boundary  layer  which  reaches  far  downstream,  and  although  its  maximum  tempera¬ 
ture  differences  are  small,  it  was  found  that  the  small  changes  in  the  mean  velocity 
profile  noticeably  altered  the  local  stability  characteristics  of  the  flow.  This  was  also 
verified  with  linear  spatial  stability  calculations  using  the  base  profiles  obtained  from 
the  Navier-Stokes  solutions. 

A  linear  transfer  function  technique  based  on  the  vorticity  response  at  the  wall 
surface  and  the  temperature  input  at  the  heater  strips  was  devebped.  The  heating 
mechanism  that  creates  the  TS  wave  packets  over  the  heater  strip  is  nonlinear  due  to 
the  inseparability  of  the  momentum  and  energy  equation.  However,  for  low  ampli¬ 
tude  forcing  the  vorticity  response  downstream  from  the  heater  strip,  the  nonlinear 
effects  can  be  assumed  small.  The  control  model  based  on  the  transfer  function 
assumes  prior  knowledge  of  the  disturbances  in  the  flow  field.  For  the  study  of  the 
applicability  of  wave  packet  attenuation  in  the  boimdary  layer,  the  uncontrolled  dis¬ 
turbance  field  was  first  determined  with  an  independent  computation.  Then,  by 
analyzing  the  wave  packet  disturbances  for  a  computation  with  no  control  applied,  a 
direct  comparison  can  be  made  with  computations,  when  control  is  applied. 

A  study  on  varying  the  heater  strip  width  revealed  that  for  a  width  of  typically 
one  half  the  wavelength  of  the  dominant  unstable  frequency  component,  the  vorticity 
response  was  maximized.  A  study  on  the  effect  of  varying  the  distance  of  the  sensor 
location  to  the  actuator  location  showed  that  the  cancellation  of  the  disturbances  was 
optimized  for  a  distance  of  approximately  two  disturbance  wavelength  of  the  main 
frequency  component  in  the  wave  packet.  For  a  distance  closer  to  the  disturbance 
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actuator  the  presence  of  the  above  mentioned  non-fluctuating  heat  spot  adversdy 
affects  the  control  efficiency.  Also,  for  distances  further  from  the  actuator  the  transfer 
function  loses  information  on  the  higher  frequency  range  due  to  the  strong  spectral 
filtering  of  the  boundary  layer. 

Numerical  simulations  of  the  control  of  the  boimdary  layer  transition  applying 
the  transf^  function  technique  showed  that  two-dimensional  and  three-dimensional 
wave  packet  disturbances  can  be  drastically  attenuated. 

For  a  simulation  with  no  control  applied,  nonlinear  interactions  of  frequency 
components  of  the  wave  packet  first  appeared  in  the  lower  frequency  range  of  the 
three-dimensional  spanwise  mode.  The  growth  rates  of  the  lower  frequency  com¬ 
ponents  of  the  three-dimensioiud  mode  ecceeded  the  linear  growth  rates  of  the  TS 
frequency  components.  The  two-dimensional  spanwise  mode,  however,  did  not  show 
the  low  frequency  interactions  of  the  three-dimensional  mode. 

Experimental  investigations  by  Gaster  and  Grant  (1975)  also  showed  the  low  fre¬ 
quency  interaction  mechanism  in  the  subharmonic  region  for  the  three-dimensional 
spanwise  mode  of  a  wave  packet  disturbance.  This  suggests  that  a  resonance 
mechanism  for  the  three-dimensional  disturbances  promotes  the  amplification  of  a 
large  band  of  frequency  components.  However,  neither  in  the  present  numerical 
simulations  nor  in  the  experiments  by  Gaster  and  Grant  could  a  specific  higher  order 
interaction  mechanism  be  identified  that  can  be  explained  by  one  of  the  fundamental 
breakdown  theories. 

To  verify  the  role  of  the  two-dimensional,  spanwise  modes  of  the  wave  packet 
on  the  amplification  mechanism  in  the  subharmoiuc  frequency  range,  a  numerical 
simulation  with  only  the  control  of  the  two-dimensional  part  of  the  wave  packet  was 
investigated.  For  in  this  case,  no  nonlinear  amplification  of  the  lower  frequency  com- 


81 


ponents  of  the  three-dimensional  modes  was  present.  Although  the  attenuation  of 
only  the  two-dimensional  mode  leaves  the  linear  amplification  of  the  three- 
dimensional  mode  intact,  the  amplitude  of  the  two-dimensional  mode  influences  the 
nonlinear  amplification  of  the  three-dimensional  mode  considerably. 

Future  invest^tions  of  the  control  of  the  wave  packet  disturbances  should  con¬ 
centrate  on  several  points.  The  eff^  of  the  m\plification  of  the  three-dimensional 
mode  on  the  attenuated  two-dimensional  mode  should  give  further  insight  on  the 
importance  of  three-dimensional  mode  control  concepts.  For  the  investigation  of  the 
later  stages  of  the  breakdown  process  to  turbulent  flow,  nonlinear  interaction  necessi¬ 
tates  the  incorporation  of  many  more  spanwise  modes.  We  believed  that  the  transfer 
function  concept  can  be  extended  to  also  include  quadratic  wave  interactions  which 
would  become  important  for  the  later  stages  of  the  transition  from  laminar  to  tur¬ 
bulent  flow.  Also,  this  transfer  function  technique  is  applicable  to  the  control  of  ran¬ 
dom  two-  and  three-dimensional  disturbances. 
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Fig.  9:  Amplitude-  and  Phase-  spectra  of  the  wall  normal  velocity  signal  at  y/A  y  =20 
at  several  streamwise  positions  for  the  case  with  no  control  applied. 
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Fig.  11:  Comparison  of  a)  the  phase  velocity  c  ^  and  b)  the  rate  of  amplification  a ,  of 
the  Navier-Stokes  calculation  (solid  line)  with  the  linear  spatial  stability  theory 


( '+'  symbols  )  at  the  streamwise  position,  x/Ax  =  140,  (Re  835),  by  ana¬ 
lyzing  the  wall  vorticity  signal  of  the  Navier-Stokes  simulation. 
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Fig.  12:  Comparison  of  a)  the  phase  velocity  and  b)  the  rate  of  amplification  a  of  the 
Navier-Stokes  calculation  (solid  line)  with  the  linear  spatial  stability  theory  ( '+' 
S3rmbols)  at  the  streamwise  position  x/Ax  =140,  (Re  ^  ^=835),  by  analyzing  the 
wall  normal  velocity  signal  at  y/Ay  =  30  of  the  Navier-Stokes  simulation. 
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Fig.  13:  Comparison  of  a)  tne  phase  velocity  and  b)  the  rate  of  amplification  a .  of 
the  Navier-Stokes  calculation  (solid  line)  with  the  linear  spatial  stability  theory 
( 'V  symbols)  at  the  streamwise  position  x/Ax  =  140,  (Re^^  =  835),  by  analyzing 
the  vorticity  signal  at  y/Ay=20  of  the  Navier-Stokes  simulation. 
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Fig.  14:  Streamwise  amplification  of  the  wall  vorticity  for  several  frequency  compo¬ 
nents  of  the  wave  packet  disturbance  with  a)  the  presence  of  a  heated  segment 
on  the  plate  surface  (steady  wall  temperature  distribution  shown  below  figure) 
and  b)  without  the  heated  segment. 
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Fig.  15:  Comparison  of  the  streamwise  amplification  of  four  frequency  components  of 
the  wave  packet  disturbance  with  a  steady  (passive)  heater  temperature  be¬ 
tween  70  ^  x/Ax  ^90  (solid  line)  and  without  a  heater  strip  (dashed  line)  for  a) 
the  waU  vorticity,  b)  the  wall  normal  velocity  component  y/Ay=10. 
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Fig.  16:  Comparison  of  the  rate  of  amplification  -  a  .  of  the  wall  vorticity  signal  at 
x/Ax=140  for  the  Ccise  with  the  passive  heater  (at  70  ^  x/Ax  <  90)  and  without 
the  passive  heater  present. 
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Fig.  17:  Steady  streamwise  velocity  and  vorticity  profiles  at  position 
x/ AxsIOO,  772  for  the  case  without  the  passive  heater  and  the  difference 
and  0)  ^  for  the  case  with  the  heater  and  the  case  without  the  passive 

heater. 
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Fig.  19:  Schematic  representation  of  the  transfer  function  technique  relating  the 
temperature  input  at  location  to  a  wall  vorticity  response  at  x^. 
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Fig.  22:  Temperature  control  signal 


T  ,  for  the  controller  strip. 
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Fig.  25;  The  rate  of  amplification,  the  phase  velocity  and  the  ratio  of  the  spanwise  to 
stream  wise  wave  number  from  the  linear  sp-tial  stability  theory  at  two 
streamwise  positions  Rej,=  800, 1200,  for  gamma  =  0,  and  four  oblique 
modes  gamma  =  10,2030,40. 
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Fig.  26;  Instantaneous  wall  temperature  at  t/At  =  70  with  spanwise  forcing  of  the  first 
heater  strip  and  a  passive  second  heater. 
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Fig.  27:  Disturbance  signal  of  the  spanwise  vorticity  several  time  levels  for  a) 

the  two-dimensional  mode,  and  b)  the  three-dimensional  spanwise  mode. 
End  of  forcing  is  at  t/At  =  120. 
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Fig.  28:  Disturbance  signal  of  the  wall  normal  velocity  at  y/Ay  =  40  at  several 
time  levels  for,  a)  the  two-dimensional  spanwise  mode,  and  b)  the 
three-dimensional  spanwise  mode.  End  of  forcing  is  at  t/At  =  120. 


x/Ax 


200 


Heater  Strip 

Fig.  29:  Streamwise  disturbance  velocity  comjwnent  at  t/At  =,  130  for  a)  the 

two-dimensional  mode  and  b)  the  three-dimensional  spanwise  mode 
with  no  control  applied. 


Fig.  31:  Spanwise  vortidty  disturbance  at  t/At  =  1200  for  a)  the  two-dimensional 
mode  ©  and  b)  the  three-dimensional  spanwise  mode  ©  .  ,  with  no 

z  k»0  * 

control  applied. 
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Fig.  32:  Amplitude-  and  Phase-  spectra  of  the  spanwise  wall  vorticity  signal  for  a)  the 
two-dimensional  mode  o)'  .  ,  and  for  b)  the  three-dimensional  mode  w' 

z  M  z  k-1 

with  no  control  applied, at  x  /Ax  =  90,100,110,120,130 . 
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Fig.  33:  Amplitude-  and  Phase-  spectra  of  the  wall  normal  velocity  signal  at  y/ Ay=40 
for  a)  the  two-dimensional  mode  and  for  b)  the  three-dimensional  mode 
v'.  with  no  control  applied,  at  x/A  x  =  90,100,110,120,130. 
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Fig.  34:  Amplitude-  and  Phase-  spectra  of  the  spanwise  vorticity  signal  with  no 

control  applied,  for  a)  the  two-dimensional  mode  oo  at  the  wall  surface 
and  b)  the  three-dimensional  mode  a'  ^  .  at  the  wall  surface,  and  for  c)  the 

Z  kal 

three-dimensional  mode  at  y/Ay  =  10,  at  the  streamwise  positions 
x/Ax  =  1 80, 260, 340, 420 , 
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Fig.  34c:  continued 


23 


Fig.  35:  Amplitude-  and  Phase-  spectra  of  the  wall  normal  velocity  signal  at  y/Ay=40 
for  a)  the  two-dimensional  mode  v'  „  and  for  b)  the  three-dimensional  mode 

k«0 

v'  with  no  control  applied,  at  x/Ax  =180,260,340,420,500. 


124 


I  I  t  I  I  I 


J2162024280  4  8  12i6  20  24  28 

friqoMUX  iitfiez 

Fwmx  4  Ttatmx  8 

noma  4  AMD  &.  TIm  pomr  qiaetm  of  two  dniMsiiooal  mvas  at  diffoant  loeatioaa:  Agon  4, 
mt—uiiMim;  fignra  5,  bud  oa  the  modeL 


n  j 


U  16  30  24  28  320 


Fhobk  6 


12  ]8  20  24  28 


twm  7 


PToinaithMPT.PiiilmiMH  ofthaipaattaafobliqiiawaa— faae»  0.wfaMM*aiMx 
P,e  X  lt.4  X  10-^:  figon  8,  aipwaiinal  wilm  figin*  7  oompofri  bahxaioT. 


Fig.  36:  Comparison  of  the  experimental  amplitude  spectra  with  the  computed  spectra 
of  the  linear  spatial  stability  theory  from  Caster  and  Grant  (1975).  The 
measured  spectra  are  taken  at  a  wall  normal  distance  outside  the  boundary 
layer  y/5  =  1.1 . 
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Fig.  38:  The  streamwise  velocity  component  u'  in  the  horizontal  plane  at  the  distur 
bance  maximum  (color  contours)  and  the  wall  temperature  disturbance  0' 
(line  contours)  at  t/At  =  70  during  the  forcing  cycle. 
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Fig.  39:  The  wall  normal  velocity  component  v '  in  the  horizontal  plane  at  the 

disturbance  maximum  (grey-shaded  contours)  and  the  wall  temperature 
disturbance  0'  (line  contours)  at  t/At  =  70  during  the  forcing  cycle. 
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Fig.  40;  The  spanwise  velocity  component  w '  in  the  horizontal  plane  at  the 

disturbance  maximum  (grey-shaded  contours)  and  the  wall  temperature 
disturbance  6'  (line  contours)  at  t/At  =  70  during  the  forcing  cycle. 
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Fig.  42:  The  wall  normal  vorticity  o)'  in  the  horizontal  plane  at  the  disturbance 

maximum  (grey-shaded  contours)  and  the  wall  temperature  disturbance  0' 
(line  contours)  at  t/At  =  70  during  the  forcing  cycle. 
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Fig.  44:  The  streamwise  vorticity  ©'  in  the  horizontal  plane  at  the  plate  surface 

(grey-shaded  contours)  and  the  spanwise  vorticity  to'  (line  contours)  at  t/At 

=  1200. 


Fig.  45:  Modulus  and  argument  of  the  transfer  function  at  several  streamwise 

positions  with  the  heater  strip  width  130  S  x/Ax^  145  fora)  the  two-dimensional 
mode  (k=0)  and  b)  the  threenJimensional  span  wise  horizontal  mode  (k=l ). 
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Fig.  46:  Amplitude-  and  Phase-  spectra  of  the  spanwise  vorticity  signal  (o  at  the 
wall  surface  for  the  two-  dimensional  mode  at  several  streamwise  positions 
with  control  applied.  The  control  strip  is  at  130  <  x/Ax  <  145. 
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Fig.  47:  Amplitude-  and  Phase-  spectra  of  the  spanwise  vorticity  signal  o)  at 
y/Ay=10  for  the  three-dimensional  mode  at  several  streamvydse  positions 
with  control  applied.  The  control  strip  is  at  130  S  x/Ax  ^145 
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Fig.  48:  Amplitude-  and  Phase-  spectra  of  the  wall  normal  velocity  signal  v'  at 
y/  Ay=40,  for  the  two-dimensional  spanwise  mode  at  several  stream  wise 
positions  with  control  applied.  The  control  strip  is  at  130<  x/Ax  <145. 


Fig.  49:  Amplitude-  and  Phase-  spectra  of  the  wall  normal  velocity  signal  v'  at 
y/Ay=40,  for  the  three-dimensional  spanwise  mode  at  several  streamwise 


positions  with  control  applied.  The  control  strip  is  at  130<  x/Ax  <145. 
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Fig.  50;  The  spanwise  vorticity  (o'  in  the  horizontal  plane  at  the  plate  surface  for  the 
control  of  the  two-dimensional  spanwise  mode  only  (grey-shaded  contours) 
and  the  spanwise  vorticity  0)'  for  no  control  (line  contours)  at  t/At  =  1200. 
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Appendix  A:  Fourier  Transform  of  Non-Periodic  Signals  and 
Estimation  of  the  Transfer  Function 

For  the  practical  implementation  of  the  discrete  Fourier  Transform  (DFT)  we  are 
restricted  to  finite  length  sequences  of  digital  signals  at  selected  values  of  the  angular 
frequency.  A  periodic  signal  of  period  N  has  the  property  that 

u(.n+N)  =  u(n).  (A.l) 

This  periodic  sequence  can  be  represented  by  the  discrete  Fourier  Transform  of  a 

finite-length  sequence  u  (0),u  (1 (N-1)  by 

u(n)=  s'  U(m )e-<2«'N)nw  ^  q  I  2,  •  •  • ,  N-\)  (A.2) 

m  =  0 

and  the  inverse  discrete  Foiuier  Transform  (IDFT)  of  a  finite-length  frequency 
sequence  U(0),U(1), . ,U(N-1)  by 

U(m)  =  ^^£u(n (m  =  0,1,2,  •  •  • ,  N-1)  (A.3) 

The  equations  A.2  and  A3  form  a  DFT  pair.  The  digital  signal  u(n)  and  the  frequency 
signal  U(m)  are  assumed  to  be  periodic  sequences,  since  they  are  both  associated 
with  Fourier  series.  Hence  if  u(n)  is  actually  zero  for  n<0  or  n>N-l  ,  the  DFT  will 
not  be  aware  of  this,  since  it  assumes  the  digital  signal  is  periodic  and  exists  for  all 
integers  n.  Thus  one  mxist  be  careful  in  interpreting  the  DFT  and  IDFT  outside  the 
range  of  0^  ^-1  and  0^  ^  -1. 

In  A.2  and  A.3  the  quantity  u(n),  n  =  0,13,  •••  N-1  is  the  real,  sampled  time-data 
function  of  the  computation  at  some  position  (xo,yo)  in  the  flow  field.  The  above  for¬ 
mula  assume,  that  the  length  of  the  sampled  data  to  be  transformed  consists  of  N 

2k 

points  spaced  At  apart.  In  our  computations,  the  timestep  is  chosen  as  At  =  — , 
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where  L  is  an  integer  number  for  proper  temporal  resolution,  and  P  is  a  typical  fre¬ 
quency  of  the  most  amplified  wave  component  in  the  wave  packet.  The  DFT  is 
defined  at  N  discrete  values  of  the  frequency  extending  up  to  the  frequency 
/s  =  NA/,  where  A/  =  l/T^  is  the  elementary  bandwidth,  and  Tp  is  the  time  period 
of  the  data  sample  (T  =  NAf).  For  a  time  waveform  to  be  faithfully  reconstructed,  it 
is  necessary  that  the  highest  frequency  present  in  the  signal  is  less  than  or  equal  to 
the  folding  frequency  /j^  =  /j/Z  In  our  computations  the  folding  frequency  is  deter¬ 
mined  by  ff4  =  pL/4jc. 

Before  the  digital  time  waveform  is  analyzed  a  data  windowing  or  tapering  is 
applied  to  avoid  a  discontinuity  at  the  end  of  the  data  set  and  to  increase  the  number 
of  the  sample  points  to  N^.  After  the  mean  flow  values  were  subtracted  from  the 
time  traces  a  linear  taper  to  the  first  and  last  10%  of  the  data  is  multiplied  by 


Pn  = 


l-cos(^) 

Ni 


l-cos(^(  n-N)) 

Ni 


where  Ni  =  N/10. 


0  ^  n  ^  Ni 
elsewhere 
N-N-i  <  n  <N 


(A.4) 


After  this  tapering,  zeros  are  added  to  the  data  samples  so  that  a  total  sample 
length  consists  of  2^^,  or  4096  elements. 

The  Fourier  transform  of  the  discrete  time  data  is  then  found  by  using  an  FFT 
algorithm  based  on  the  DFT  pair  A.2  and  A.3.  However  one  problem  always  arises 
when  estimating  the  discrete  spectra  because  of  the  finite  length  of  the  data.  To 
represent  the  finite  data  length  mathematically,  the  actual  computationally  generated 
data  sequence  Ucomp  can  be  expressed  in  terms  of  a  sequence  of  infinite  data  length 
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ujin )  multiplied  by  a  window  function  g{n)  so  that 

“comp(« )  =  )g(n)  (A.5) 

with 

fl  «=0,1, 

g  (n )  -  |q  elsewhere  (A.6) 

In  the  Fourier  transformed  space  the  equation  A.5  is  expressed  by 

L/comp^"*)  =  *  G(m)  m=0,l,  •  ■  •  >1-1  (A.7) 

where  U»  denotes  the  Fomier  transform  of  the  infinite  data  length  record,  m  is  the 
frequency  index,  and  the  asterisk  denotes  the  convolution.  The  equation  A.7  indi¬ 
cates  that  the  computational  Fourier  transform  is  given  by  the  true  Fourier  transform 
and  its  convolution  with  the  window  function  in  the  frequency  space.  From  the  pro¬ 
perties  of  the  convolution  integral,  it  is  known  that  Ua^  will  approximate  U.,  only 
in  the  case  when  G(m)  is  confined  to  a  narrow  band  of  frequencies.  In  the  limit, 
where  G(m)  approximates  a  delta  function,  Ua^  =  U„.  The  wave  packets  in  our 
computation  consist  of  a  nonperiodic  signal  which  disappears  once  the  disturbance 
has  completely  passed  the  recording  station,  thus  leakage  of  frequency  components 
through  the  sidebands  of  the  window  function  (the  taper  function)  plays  a  minor 
role.  Instead,  pruning  of  the  nonperiodic  signal  by  adding  zeros  to  the  time  sequence 
greatly  sharpens  the  peak  of  the  box  car  function  A.6  and  improves  the  spectral  esti¬ 
mate  considerably.  Finally,  the  amplitude  spectrum  A  {m )  and  the  phase  spectrum 
(|>(m)  of  the  Fourier  transformed  time  data  are  found  through  the  relationships 

A(m)  =  ,  pirn)  =  tan"^(-^).  (A.8) 

where  and  b„  are  the  real  and  imaginary  parts  of  U  (m ),  respectively. 
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The  general  relation  between  input  and  output  in  continuous  time  of  a  linear 
system  is  given  by 

m 

y(f)  =  J  k(x)x{t-x)  dx  (A.9) 

For  our  computations  we  restrict  the  above  integral  to  a  physically  realizable 
system  by  assuming  that  the  input  to  the  system  for  x<0  is  zero.  Then,  by  taking  the 
Fourier  transform  on  both  sides  of  A.9,  one  arrives  at: 

Y(m)  =  H(m)X(m)  (A.IO) 

where  Y(m)  and  X(m)  are  the  discrete  Fourier  Transform  of  the  input  and  output  sig- 

ncils  respectively. 

From  the  relation  A.IO,  one  can  easily  find  the  transfer  function  H(m).  The  rela¬ 
tion  A.10  is  strictly  true  only  in  the  absence  of  any  noise.  In  our  numerical  data, 
there  will  always  be  some  computational  noise  present.  This  noise  level  is  deter¬ 
mined  by  the  accuracy  of  the  machine  and  the  convergence  criteria  used  in  our 
numerical  scheme.  The  convergence  criteria  is  set  so  that  the  magnitude  of  the  trun¬ 
cation  is  several  orders  lower  compared  to  the  level  of  flow  perturbation. 

The  estimation  of  H(m)  in  A.14  is  limited  to  the  range  of  frequencies  where  nei¬ 
ther  X(m)  or  Y(m)  is  exactly  zero.  To  avoid  this  mathematical  dilemma  of  dividing  by 
zero,  a  filter  in  the  frequency  domain  is  applied  on  H(m)  whenever  X(m)  or  Y(m)  is 
below  the  machine  accuracy.  Below  this  threshold,  the  convolution  loses  all  informa¬ 
tion,  and  reconstructing  that  frequency  component  becomes  impxsssible.  Since  the 
transfer  function  gradually  falls  off  to  zero  for  higher  frequency  components,  a  sim¬ 
ple  step  function  filter  is  used  when  the  threshold  of  machine  accuracy  is  reached. 


Appendix  B:  Finite  Difference  Approximations  for  the 

Energy  Equation 
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The  integration  domain  is  divided  into  N  grid  points  in  the  x-direction  and  M 
grid  points  in  the  y-direction.  The  grid  is  uniform  in  both  directions.  In  the  following 
expressions,  the  subscripts  m  and  n  refer  to  constant  grid  lines  in  the  y-  and  x- 
direction  respectively.  The  subscript  p  denotes  the  time.  As  a  result,  the  spatial  coor¬ 
dinates  arex  =  nAz  and  y  =  mAy  and  the  time  ist=  pAt. 

The  first  derivative  in  time  for  p  >  1  at  and  is  approximated  by: 

do  1 

~  "2^^  ^  ~  4  "*■  0(At^)  (B.l) 

and  for  p=l: 

~  ~^(^H^,p  “  +0(Af)  (B.2) 

where  6  is  the  temperature. 

For  the  following  spatial  derivatives,  the  superscript  p  for  the  time  has  been  dropped, 
and  each  derivative  is  evaluated  at  the  current  time  step.  Also  the  coordinate  z  refers 
to  the  either  the  streamwise  (x)  or  the  wall  normal  direction  (y),  resp)ectively. 

For  2^n  ^  -2  and  2^  ^  -2: 

1 

0k-2  -  8  0k-i  +  8  01^1  -  0  )  +  O  (Az")  (B.3) 

^'k  =  -  ek-2  +  16  0k_i  -  30  0*  +  16  0t^i  -  0k^2  )  +  O  (Az^)  (B.4) 

at  m=1  (the  plate  surface),  the  y-derivatives  are  given  by: 

■^lt  =  :j^(-AF-650o  +  40, +  6602  -  403  -  04  )  +  O  (Ay^) 


(B.5) 
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■^li  = - +  145  00  -  304  01  +  174  0,  -  16  0,  +  0.  )  +  O  (Ay^)(B.6) 

120Ay2  u  i  ^  3  4 

AF  is  a  known  scalar  at  the  wall  surface  at  the  timestep  p: 


AF  =  12Ay2  •  Pr 


d0  , 
dt 


1  a^0 

Re  Pr  dx^ 


'waU 


(3.7) 


At  the  outflow  boundary,  at  grid  lines  N  and  N-1,  the  first  and  second  derivatives 
are  approximated  by: 


-  4  0N_3  +  27  0|^_2  -  108  0,^_i  +  85 

-  )  +  0(Ar*)  (B.8) 

do  1 

=  50^^  ~  ^  ®W-4  +  60  0^/_3  -  120  0>/_2  +  65  0>/_i 

+  12  0^  )  +  O(Ax^)  (B.9) 


^  ^  ^  ^  +  26  0/^  (B.IO) 

+  Zk^AxHBn-Bn^^  )  +  CXAr^) 

At  the  grid  line  next  to  the  inflow  boimdary  n=l,  and  the  grid  line  next  to  the  upper 
botmdary  m  =  M-1,  the  following  one-sided  approximations  are  used  for  the  y- 
derivatives 


d0  1 

-  ‘|22iy  ^  ^  ®M+10  0M-1~18  0M-2+6  0m -3“  8m-4  )  +  O  (Ay^)  (B.ll) 
eP’Q  1 

^  1O8m~150m-i-40m-2+140m-3-60m-4+Om-5  )  +  O  (Ay^)  (B.12) 

The  temperature  dependent  terms  of  the  vorticity  transport  equation  (2.3)  are  refor¬ 
mulated  with 


Ri  = 


ax  Bx 


(B.13) 
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p  -  ^  -  w  ^ 
/•  >2 


R3=^=»" 

^  dx^ 


V'" 


ae 

ax 

30 


,  a^e 

■"'"'ax^ 


12 


+V 


/  a^6 


a^ 


ae  30  . . .  a^e 


Ks=  — ^  =  t>  +1> 


(B.14) 


(B.15) 


(B.16) 


(B.17) 


ay2  “  ax  ^  axay 
The  prime  on  the  viscosity  v  indicates  a  derivative  with  respect  to  temperature  0, 
which  is  known  analytically  from  equation  (Z7). 

With  the  above  approximations  the  temperature  dependent  terms  can  be  evaluated: 


^  ‘  0n-l,m+8  )  (B.18) 

0n,m-2~8  0n^-l+8  0n,m-l-  ^H/n-2  ) 

14^2{®«-2*»«“®  0fi-l^+8  0n+l^-  %+2^^  (B.21) 

To7~2  ^  *”  0„^+16  0,1+1^—  07i+2,m  ) 

izAx 

^4^»i,»n  ~  ®n^+l~  0n^+2j  (B.22) 

.  2  ^  ®ii^-l~30  0„^+16  0n^+l“  0n,m+2  ^ 

12Ay‘ 

^5^n,m  ~  0n^+l“  0n,m+2)  (B.23) 

^®n-2,m~0  0»i+l,in“  0n+2^) 

lii/icAy  ^  0n-l^-2+8  0„+l^-2-  0«+2^-2 

“®fi-2,«-l+^  0>i+l.m-l+®  ®>i+2,m-2 

+®  0«-2,m+l~04  0„_i^+l+64  0„+i^+l-8  0„+2^+l 
+  0n-2,M+2+®  ®n-l,jn+2“®  0n+l,fn+2+  0n+2,m+2  ^ 
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